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ABSTRACT 


Measurements  were  made  of  the  heat  loss  and  recovery  tem¬ 
perature  of  a  fine  hot-wire  at  a  nominal  Mach  number  of  5.  8.  Data 
were  obtained  over  an  eight-fold  range  of  Reynolds  numbers  in  the 
transitional  regime  between  continuum  and  fr ee-molecule  flow.  At 
high  Reynolds  numbers,  the  heat  transfer  data  agree  well  with  the 
results  of  JLaufer  and  McClellan,  which  were  obtained  at  lower  Mach 
numbers.  At  lower  Reynolds  numbers,  the  results  indicate  a  monotonic 
transition  between  continuum  and  free  molecule  heat  transfer  laws. 

The  slope  of  the  heat  transfer  correlation  also  appears  to  vary  mono- 
tonically,  with  Nu  / Re  at  high  Reynolds  numbers  and  Nu  Re  for 
Re  <  <  1. 

Data  on  the  wire  recovery  temperature  (corresponding  to  zero 
net  heat  transfer)  were  obtained  for  free- stream  Knudsen  numbers 
between  0.  4  and  3.  0.  Comparison  with  previous  data  suggests  that  for 
Mach  numbers  greater  than  about  two  the  normalized  variation  of 
recovery  temperature  in  the  transitional  regime  is  a  unique  function 
of  the  free- stream  Knudsen  number. 

The  steady- state  hot-wire  may  be  used  to  obtain  two  thermody¬ 
namic  measurements:  the  rate  of  heat  transfer  from  the  wire  and  the 
wire  recovery  temperature.  An  illustrative  experiment  was  performed 
in  the  wake  of  a  transverse  cylinder,  using  both  hot-wire  and  pressure 
instruments  in  a  redundant  system  of  measurements.  It  was  shown  that 
good  accuracy  may  be  obtained  with  a  hot-wire  even  when  the  Reynolds 
number  based  on  wire  diameter  is  small. 
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I.  INTRODUCTION 

Hot-wire  probes  have  been  used  for  a  number  of  years  to 
measure  the  mean  properties  of  a  fluid  stream.  The  hot-wire  is 
capable  of  providing  two  thermodynamic  measurements:  the  rate  at 
which  heat  is  transferred  from  the  wire  to  the  stream  and  the  wire 
recovery  temperature  for  zero  heat  transfer.  If  sufficiently  accurate 
relations  are  known  for  the  heat  loss  and  recovery  temperature  as  a 
function  of  the  properties  of  the  flow,  only  one  additional  thermodynamic 
measurement  is  necessary  to  specify  the  flow  field  uniquely. 

The  present  investigation  was  undertaken  with  two  purposes  in 
mind.  First,  in  using  the  hot-wire  to  measure  mean  quantities  in  low 
density  compressible  flow,  it  is  necessary  to  have  an  accurate  calibration 
of  the  wire  heat  loss  and  recovery  temperature  as  a  function  of  Mach 
number  and  Reynolds  number.  For  high  Reynolds  number  flows,  these 
relations  are  well  represented  by  the  data  of  Laufer  and  McClellan*.* 
Often,  however,  it  is  necessary  to  use  the  hot-wire  in  regions  where 
the  Reynolds  number  based  on  wire  diameter  is  small.  Previous 
investigations^'  ^  have  illustrated  the  general  qualitative  features 
of  the  Nu  -  Re  relation  at  small  Reynolds  numbers,  but  the  experimental 
uncertainties  are  too  large  for  calibration  purposes.  The  first  objective, 
therefore,  was  to  obtain  an  accurate  hot-wire  calibration  at  high  Mach 
number  in  the  transition  regime  between  continuum  and  free  molecule 
flow. 


*  Superscripts  denote  references  at  the  end  of  the  text. 
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The  second  purpose  of  this  investigation  was  to  provide  further 

information  on  two  features  of  low  Reynolds  number  flow  which  are 

important  from  a  theoretical  point  of  view.  One  of  these  is  the  recovery 

temperature  phenomenon  which  was  theoretically  and  experimentally 

8 

investigated  by  Stalder,  Goodwin,  and  Creager  „  It  was  found  by  these 
authors  that  the  equilibrium,  or  recovery, temperature  of  an  unheated 
cylinder  in  free-molecule  flow  is  greater  than  the  stagnation  temperature 
of  the  stream.  Since  the  high  Reynolds  number  equilibrium  temperature 
is  less  than  stagnation,  the  recovery  ratio  (Taw/TQ)  serves  as  one 
measure  of  the  transition  from  continuum  to  free  molecule  flow.  In 
the  present  investigation,  measurements  were  made  of  the  wire  recovery 

temperature  in  the  range  of  Knudsen  numbers  between  0.  4  and  3.  0. 

1  2 

These  and  previous  ‘  data  provide  an  accurate  high  Mach  number 
relation  between  the  recovery  ratio  (Taw/TQ)  and  Kn  over  the  complete 
range  from  continuum  to  free  molecule  flow. 

At  high  Reynolds  numbers,  the  heat  transfer  to  the  stagnation 
point  of  a  blunt  body  may  be  accurately  predicted  using  boundary  layer 
theory.  Recently,  a  considerable  amount  of  attention  has  been  given  to 
the  problem  of  determining  appropriate  corrections  to  the  high  Reynolds 
number,  high  Mach  number  boundary  layer  solutions  which  would  extend 
their  validity  to  lower  Reynolds  numbers.  Several  alternative  formulations 
of  this  problem  have  been  given  in  the  literature,  and  more  experimental 
evidence  is  needed  to  evaluate  the  results.  The  data  obtained  in  this 
investigation,  when  combined  with  the  results  of  Laufer  and  McClellan, 
show  the  departure  of  the  Nu-Re  relation  from  the  high  Reynolds 
number  correlation  and  the  approach  to  free-molecule  flow  for  an 
infinite,  conducting  cylinder. 
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In  many  physical  situations,  it  is  difficult  to  obtain  precise 
data  using  standard  pitot  probes.  Particularly  in  regions  of  low 
density,  instrumentation  sufficiently  large  to  obtain  reliable  results 
may  introduce  appreciable  disturbances  into  the  flow.  Total  temperature, 
mass  flow,  and  velocity  profiles  in  the  wake  of  a  transverse  cylinder 
are  reported  here  to  illustrate  the  use  of  the  hot-wire  in  such  regions. 

By  using  both  hot-wire  and  pressure  measurements,  the  flow  field 
was  over-determined,  and  this  redundancy  permitted  an  estimate  of  the 
accuracy  of  the  hot-wire  measurements. 
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II.  EXPERIMENTAL  METHOD 

II.  1.  Equipment 

The  measurements  were  performed  in  Leg  1  of  the  GALCIT 
hypersonic  wind  tunnel  at  a  nominal  Mach  number  of  5.  8.  This  tunnel 
has  a  3  in.  by  3  in.  core  of  uniform  flow.  Other  features  of  this  facility 
are  described  in  Reference  9.  The  stagnation  pressure  was  varied 
between  0  and  95  psig.  at  a  nominal  stagnation  temperature  of  226°F. 
Because  of  the  tunnel  wall  boundary  layer,  the  actual  test  section  Mach 
number  varied  with  stagnation  pressure  as  shown  in  Figure  1.  Figure  2 
shows  the  variation  of  stagnation  temperature  with  pressure  which  was 
introduced  by  the  automatic  temperature  regulation  system. 

Under  certain  conditions,  free- stream  turbulence  may  effect 
the  measured  heat  transfer  rate,  Demetriades^  found  that  the  free- 
stream  integrated  r.  m.  s.  mass-flux  and  temperature  fluctuations  in 
Leg  1  were  on  the  order  of  0.  4  per  cent.  The  dominant  turbulence 
scale,  given  by  the  velocity  divided  by  the  frequency,  was  on  the  order 

3  5.. 

of  10  to  10  times  the  wire  diameter,  and  by  the  criteria  of  Reference 
12  should  have  caused  no  measurable  effect  on  the  heat  loss.  Additional 
quantitative  turbulence  measurements  are  necessary  before  a  similar 
statement  can  be  made  about  the  wake  results. 

To  determine  the  local  Mach  number,  pitot  pressure  measurements 
were  made  within  0.  5  in.  of  the  hot-wire  with  the  .  063"  I.  D.  stainless 
steel  probe  shown  in  Figure  3,  and  the  local  Mach  number  computed 
assuming  isentropic  flow  between  the  reservoir  and  test  section.  The 
stagnation  pressure  was  measured  by  a  total  head  tube  directly  upstream 
of  the  throat.  Traverses  normal  to  the  flow  direction  showed  that  the 
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local  Mach  number  at  the  wire  was  identical  to  that  measured  by  the 
probe  within  0.  5  per  cent. 

Stagnation  temperatures  were  measured  1  inch  ahead  of  the 
nozzle  throat  with  a  Pratt  and  Whitney  type  iron-constantan  thermocouple, 
and  all  computations  were  made  assuming  isentropic,  adiabatic  expansion 
to  the  measured  test  Mach  number. 

The  detachable  hot-wire  holders  consisted  of  a  conical  brass  body, 
a  thin  double  wedge  brass  support  blade,  and  two  sewing  needles  to  hold 
the  wire.  The  needles  were  attached  to  and  electrically  insulated  from 
the  brass  blade  by  epoxy  resin;  Figure  3a  illustrates  the  holder  con¬ 
figuration.  The  complete  hot-wire  assembly,  including  holder,  sting 
support,  and  total  pressure  tube  is  shown  in  Figure  3b.  The  hot-wire 
was  soft- soldered  to  the  needle  supports  within  .  002"  of  the  tips; 
a  typical  installation  is  shown  in  Figure  4.  Details  of  the  electrical 
connectors  and  head  assembly  may  be  found  in  Reference  10. 

Figure  5  is  a  drawing  of  the  calibration  oven  used  to  determine  the 
temperature-resistance  relation  of  the  wires.  Heat  was  supplied  by  a 
helically  wound  nichrome  ribbon,  and  sufficient  thermal  lagging  was 
used  to  make  the  temperature  gradients  within  the  oven  negligible.  An 
iron-constantan  thermocouple  placed  0.  25  inch  from  the  wire  gave  the 
oven  temperature.  Measurements  were  made  only  when  the  oven  had 
reached  thermal  equilibrium  (temperature  change  less  than  0.  2°F  per 
minute). 

The  electrical  instrumentation  consisted  of  a  precision  Leeds 
and  Northrup  "K-2"  potentiometer  accurate  to  one  part  in  5,  000,  five-45 
volt  dry  cells  to  supply  the  wire  current,  a  variable  series  resistance 
for  current  control,  and  an  auxilliary  milliammeter  and  bridge  as  shown 
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in  Figure  6.  Measurements  of  the  potential  drop  across  the  wire  and 
the  standard  one  ohm  resistor  were  made  with  the  K-2  potentiometer, 
giving  the  wire  current  and  resistance.  All  readings  were  corrected 

V 

for  the  small  line  resistance  between  the  hot-wire  and  the  potentiometer. 
During  each  series  of  measurements,  the  current  was  constant  within 
0.  02  per  cent. 

II.  2.  Measurements 

To  determine  the  Nusselt  number  and  recovery  temperature  of 

the  wire,  the  voltage  drop  across  the  wire  and  across  a  standard  1  ohm 

resistor  were  measured  for  several  values  of  wire  current.  The 

2 

resistance  was  then  plotted  against  i  to  determine  the  slope  and 

intercept  at  zero  current.  Extreme  care  had  to  be  exercised  in 

achieving  constant  tunnel  stagnation  temperature,  since  a  change  of  1°F 

2 

was  easily  detectable  in  the  R  vs.  i  plot. 

In  performing  the  wake  test,  a  continuous  curve  of  wire  temperature 
vs.  position  in  the  wake  is  desired  for  a  fixed  wire  current,  and  a 
Mosely  Autograph  x-y  plotter  was  used  for  this  purpose.  The  traverse 
mechanism  was  geared  to  a  linear  resistance  coil,  and  a  small  constant 
current  through  this  resistance  created  a  voltage  drop  proportional  to 
the  traverse  position.  This  signal  and  the  voltage  drop  across  the  hot¬ 
wire  were  the  two  electrical  inputs  for  the  plotter.  Scale  factors  were 
obtained  by  measururing  the  voltage  drop  across  the  wire  with  the  K-2 
potentiometer  at  two  known  values  of  (y/D). 
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II.  3.  Wire  Calibration  and  Annealing 

The  resistance  of  each  wire  was  measured  at  several  temperatures 
in  the  calibration  oven.  At  the  end  of  each  series  of  runs,  the  wire 
resistance  at  room  temperature  was  again  measured;  if  the  resistance 
had  changed  more  than  a  few  tenths  of  a  per  cent,  the  data  were  dis¬ 
carded.  Several  wires  were  calibrated  both  before  and  after  prolonged 
testing  to  determine  the  effect  of  annealing  and  exposure  on  the 
temperature- re sistivity  coefficient. 

All  data  were  obtained  using  0.00010  inch  diameter  Pt  -  10  per 
cent  Rh  wire  from  a  single  spool.  (Some  preliminary  tests  using 
0.00005  inch  diameter  wire  are  discussed  in  Appendix  D.  )  No  attempt 
was  made  to  measure  the  wire  diameter  directly,  since  it  was  felt  that  the 
manufacturer's*  drawing  procedures  were  of  greater  accuracy  than  any 
available  optical  measurement  technique.  However,  each  wire  was 
examined  visually  under  a  microscope  to  detect  any  macroscopic 
irregularities.  The  mounted  wire  resistance  and  length  were  also 
carefully  measured  and  checked  against  the  manufacturers  quoted 
resistance  of  11,  580  ohms/ft.  at  68°F.  The  mean  deviation  of  the 
unannealed,  installed  resistance  from  this  reference  value  was  less 
than  1  per  cent,  indicating  high  uniformity  of  the  wire. 

Through  prolonged  use,  the  wire  resistance  increased  between  2 
and  8  per  cent.  Since  it  was  necessary  to  limit  the  increase  between 
successive  calibrations  to  a  few  tenths  of  a  per  cent,  several  special 
tests  were  conducted  to  determine  the  major  factors  involved.  It  was 
found  that  a  large  part  of  the  resistance  shift  could  be  eliminated  by 


*  Sigmund  Cohn,  Mt„  Vernon,  New  York. 
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heating  the  wire  to  a  dull  glow  for  several  minutes  prior  to  tunnel  use. 

This  result  is  in  substantial  agreement  with  the  experience  of 

Spangenberg  (c.  f.  ,  Reference  18,  Figure  9).  To  check  the  influence  of 

0 

starting  loads,  one  wire  (  =  325)  was  subjected  to  seven  consecutive 

tunnel  starting  cycles  at  20  psig  while  immersed  deep  within  the  sidewall 
boundary  layer.  The  resulting  resistance  change  was  less  than  0.  1  per 
cent. 

It  can  be  inferred  that  the  primary  causes  of  resistance  change 
are  associated  with  exposure  to  the  free- stream  pressure  loads  while 
the  wire  is  at  an  elevated  temperature.  Further  observation  showed 
that  the  two  primary  factors  were:  (1)  large  step-wise  increases  in 
resistance  if  the  previous  maximum  temperature  and  pressure  of 
operation  were  exceeded  [^"deformation"  phenomenon]  ;  (2)  a  much 
smaller  increase  which  was  a  function  of  time  Jj'creep"  phenomenon]  . 
Since  all  wires  were  initially  normal  to  the  flow  direction,  a  certain 
amount  of  stretching  was  to  be  expected.  Figure  4  shows  a  typical 
deformation  caused  by  an  overheat  of  about  0.  8  at  90  psig  stagnation 
pressure.  Higher  overheats  resulted  in  larger  wire  curvature  and  were 
not  used  in  this  investigation.  The  creep  phenomenon  was  negligible 
after  the  first  few  hours  of  tunnel  exposure. 

Two  factors  limited  the  aspect  ratio  of  the  wires.  First,  it  was 
found  that  wire  curvature  increased  with  aspect  ratio.  Second,  for  aspect 
ratios  greater  than  375,  wires  of  0.  00010"  diameter  were  easily  broken 
during  the  tunnel  starting  process.  Since  it  was  desired  to  keep  the  end 

<  i  < 

loss  corrections  as  low  as  possible,  the  range  320  =  =  360  was 

chosen  as  an  acceptable  compromise. 

The  temperature  coefficient  of  resistivity  ar  varied  between  0.  88 
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and  0.  93  (0F)"*  at  0°F,  with  a  mean  value  of  0.  91 .  Equilibrium  wire 
temperatures  were  less  than  270°F,  and  resistance  varies  linearly  with 
temperature  in  this  range.  Calibration  data  were  taken  between  72 
and  200°F,  and  the  measured  resistivity  coefficient  was  accurate 
to  about  1.  5  per  cent  for  each  wire.  Several  tests  indicated  that  the 
resistivity  coefficient  was  not  affected  by  tunnel  exposure,  provided  that 
the  wire  had  been  initially  annealed  to  a  dull  glow  in  still  air. 

Under  certain  conditions,  reversible  strain  gage  effects  may 
introduce  systematic  errors  into  the  equilibrium  temperature  measure¬ 
ments.  Using  the  relations  of  Morkovin  (Reference  11,  p.  24),  the 

resistance  error  /S  for  this  experiment  caused  by  strain  gauge  effects 
.  _4  _5 

would  be  (ARw/Rw)  ~  10  to  10  ,  or  considerably  less  than  the 

uncertainty  of  the  wire  calibration.  A  more  serious  error  could  be 
introduced  by  differential  thermal  expansion  between  the  wire  and  the 
support  head.  However,  no  evidence  of  this  effect  was  present  in  the 
temperature  calibrations,  and  since  tunnel  operation  was  less  severe 
from  this  standpoint  than  the  oven  calibration,  it  was  assumed  that 
thermal  strain  was  negligible. 

In  Appendix  A,  the  correction  factors  for  heat  loss  to  the  wire 

supports  are  derived.  This  analysis  is  carried  through  with  the  support 

temperature  as  a  free  parameter>and  shows  that  for  large  end  loss 

correction  the  support  temperature  has  an  important  effect  on  both  the 

Nusselt  number  and  recovery  factor.  Depending  upon  the  tip  shape, 

thermal  conduction,  and  free- stream  Reynolds  number  the  effective 

value  of  Y\  =  (Tg/TQ)  should  lie  between  .  85  and  1.  0,  and  should  vary 

with  the  wire  attachment  location.  Since  the  recovery  temperature 

correction  was  quite  sensitive  to  Y]  ,  a  special  test  (see  Appendix  B) 

s 
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was  conducted  to  measure  this  parameter.  The  value  of  /"j  was  found 
to  be  constant  over  the  complete  range  of  tunnel  conditions  encountered, 
and  equal  to  0.  903  -  .  005.  This  value  is  approximately  equal  to  the 
calculated  laminar  boundary  layer  recovery  factor  for  flow  along  the 
needle  downstream  of  the  stagnation  region. 


11 


III.  DATA  REDUCTION 

In  previous  hot-wire  investigations,  many  attempts  have  been 
made  to  secure  a  universal  correlation  of  Nusselt  number  as  a  function 
of  Reynolds  number  and  Mach  number.  Since  hot-wire  data  have  been 
obtained  in  several  different  gases  over  a  Reynolds  number  range  of 
eight  orders  of  magnitude,  for  Mach  numbers  between  zero  and  six  and 
for  wire  overheats  T  between  0  and  2,  it  is  by  no  means  a  trivial 
matter  to  find  correlations  which  represent  the  data  over  any  large 
range  of  parameters. 

The  largest  portion  of  heat  transferred  from  a  blunt  body  is 
in  the  vicinity  of  the  forward  stagnation  point.  At  high  Reynolds  numbers 
and  Mach  numbers,  the  appropriate  temperature  for  evaluating  viscosity 
and  thermal  conductivity  is  that  existing  behind  a  normal  shock  at  the 
free  stream  Mach  number.  At  low  Reynolds  numbers,  a  distinct  shock 
wave  does  not  appear;  under  these  conditions,  the  results  of  free  molecule 

g 

theory  may  be  conveniently  expressed  using  fluid  properties  evaluated  at 
the  stagnation  temperature.  At  high  Mach  numbers,  the  temperature 
behind  a  normal  shock  is  close  to  the  stagnation  temperature  and  the 
two  evaluations  will  give  equivalent,  but  not  identical,  correlations. 

In  this  paper,  the  heat  transfer  data  are  correlated  using  the  Nusselt 
and  Reynolds  numbers 

NUo  =  <hd/ko)  ,  Rec  =  (p„U„d)/A  ,  (1, 

where  /U. q  and  kQ  are  the  air  viscosity  and  thermal  conductivity 
evaluated  at  stagnation  temperature.  This  form  offers  distinct 
advantages  when  the  hot-wire  is  used  as  a  measuring  instrument  in 
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non-uniform  flow  fields.  * 

All  hot-wires  of  finite  length  lose  a  certain  amount  of  heat  to 

their  supports.  The  quantity  Nuq  represents  the  heat  transfer  rate  for 

a  wire  of  infinite  length,  and  is  related  to  the  measured  Nusselt  number 

Nu  by  the  relation 
m  3 


Nu 


=  ■f  Nu  =  5?  r _ i_ 


m 


ko  <Tw  "  Tawm^ 


] 


where  "p  _ T  is  the  Nusselt  number  end  loss  correction  factor,  i^R 

N  m 

is  the  measured  Joule  heating,  and  (T  -  T  )  is  the  measured 

'  w  awm' 

temperature  potential.  Similarly,  the  recovery  ratio  Y{  *  for  an 
infinitely  long  wire  may  be  written 


(2) 


ft  *  =  <Taw*/To>  *  ^  b  (Tawm/T„) 

5*  i. 


R  *m 


R  '  awm'  o' 


(3) 


where  ir  is  the  recovery  ratio  end  loss  correction  factor.  In 
Appendix  A,  the  quantities  and  are  derived  in  detail. 

They  depend  primarily  on  the  wire  aspect  ratio  and  measured  Nusselt 
number;  for  large  aspect  ratios  and  Nusselt  numbers,  ^  and  3? 
are  close  to  unity. 

The  present  measurements  were  conducted  using  very  small 
heating  currents.  For  i  — 0  and  a  linear  temperature  -  resistance 
relation. 


R 


•2  „ 
i  R. 


m 


=  a  R  R  ( 


(T‘  1  -  T  )  "  14 r  ^m  iXr  '  -  K 

'  w  aw'  m  < 


a  R  R 
^  r  awm  r 


awm 


(dR/di2^ 


and  Eq.  2  takes  the  simplified  form 


*  Morkovin 
more  detail. 
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and  Collis  and  Williams 
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discuss  this  problem  in 
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Nu  =  ( 


Rawm  Rr 
(dR/di^).=o 


(2a) 


2  2 
A  plot  of  R  vs.  i  (e.  g.  Figure  11)  allows  the  quantity  (dR/di  )^_Q 

to  be  determined  accurately. 

For  high  Reynolds  numbers,  the  recovery  ratio  Y\c=  <T  /T  ) 

has  been  shown  experimentally  to  be  a  function  only  of  Mach  number^; 

decreases  from  1 . 0  at  M  =  0  to  0.  95  for  Mach  numbers  greater  than 

about  two.  This  behavior  is  explained  by  the  well-known  hypersonic  "freeze" 

or  Mach  number  independence  principle, which  recognizes  that  the  viscous 

and  inviscid  flow  fields  become  independent  of  Mach  number  if  the  Mach 

number  is  sufficiently  large.  For  blunt  bodies  such  as  the  cylindrical 

wire  considered  here,  changes  in  the  flow  field  are  small  beyond  a 

Mach  number  of  three.  *  A  value  Y1  =  0.  950  was  used  as  a  reference 

c 

for  the  present  experiments  and  the  data  of  Laufer  and  McClellan*  and 
Sherman^.  Stalder  et  al^  found  =  0.  96,  and  this  value  was  used  in 

interpreting  their  results. 

The  high  Reynolds  number  relations  described  above  represent 
the  limit  of  continuum  flow.  For  free  molecular  flow,  the  heat  loss 
and  recovery  temperature  of  an  infinite,  perfectly  conducting  cylinder 
with  no  radiative  cooling  are  given  in  Reference  8: 


14 

*  Tewfik  and  Geidt  have  measured  the  angular  variation  of 
heat  transfer  coefficient  and  recovery  temperature  for  a  cylinder  at 
Mach  numbers  of  1.  3  and  6.  0.  These  distributions,  which  are  similar 
to  those  obtained  at  low  speeds,  may  be  used  to  calculate  the  temperature 
Taw  a  cylinder  of  infinite  thermal  conductivity  would  assume  for 

no  net  heat  transfer.  As  expected,  the  resulting  values  of  T  (M) 
agree  well  with  the  hot-wire  measurements. 
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Nu 


■<y-p 


2(t) 


T/T 


a  Re  Pr  ( 

o  o 


g(Sj) 


(4) 


n, 


(T  /T  )  =  (1  + 

'  aw'  o'  '  2 


)  y/z  m 


2  f(si) 

M)  w 


(5) 

(6) 


where  f(s^)  and  g(sj)  depend  only  on  the  molecular  speed  ratio 
and  the  number  of  excited  degrees  of  freedom  of  the  medium.  These 

functions  are  tabulated  in  Reference  8  for  both  monatomic  and  diatomic 

\ 

perfect  gases.  As  M  — ►-  cc  ,  Eqs.  (4)  and  (5)  approach  the  limiting  forms 


and 


( 


8'+  1 

TVr 


) 


2 

y  + 


) 


a  Re  •  Pr 
o  o 


(4a) 


(5a) 


As  in  the  case  of  continuum  flow,  the  Nusselt  number  and  recovery  ratio 
become  independent  of  Mach  number  as  the  Mach  number  increases.  * 
The  fact  that  the  recovery  ratio  in  continuum  and  free  molecule 
flow  obeys  the  Mach  number  independence  principle  suggests  that  the 
normalized  recovery  ratio 


n 


* =  7irr7r 


(7) 


is  a  unique  function  of  the  Knudsen  number  for  M  sufficiently  large.  At  lower 

Mach  numbers,  ^  would  depend  on  additional  parameters  involving  M  and  Re. 

7 

Cybulski  and  Baldwin  report  that  for  subsonic  flow  the  continuum  recovery 
ratio  persists  to  quite  high  Knudsen  numbers. 


*  Reference  15  gives  several  additional  examples  of  the  Mach 
number  independence  principle  for  both  continuum  and  free  molecule  flow. 
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All  computations  of  flow  properties  were  made  using  =  (7/5) 
and  the  Sutherland  viscosity  formula  (with  T  expressed  in  °R) 

,3/2 


M  ~  2*  270  198.  6  +  T  X  10 


lb  -  sec 

— 


(8) 


ft 


Values  of  air  thermal  conductivity  were  obtained  from  Hilsenrath  and 
Beckett^,  and  the  wire  thermal  conductivity  was  taken  from  Reference 
17,  with  the  data  linearly  extrapolated  to  the  measured  mean  wire 
temperature. 

A  derivation  of  the  end  loss  corrections  appropriate  to  transitional 
flow  is  given  in  Appendix  A.  This  analysis  assumes  that  the  wire 
temperature  is  radially  uniform,  wire  resistance  is  a  linear  function  of 
temperature,  and  the  heat  transfer  coefficient  h  =  Nuq  (kQ/d)  is  constant 
over  the  wire  length.  For  convenience,  the  supports  are  taken  to  be 
equal  in  temperature,  although  this  restriction  is  not  essential.  The 

relation  between  the  measured  Nusselt  number  Nu _  and  the  value  Nu_ 

for  an  infinite  wire  may  be  expressed  as 

(1  +  st*) 


m 


< 


t*  “  £ 


) 


*  &  (1  +  St) 

while  the  recovery  ratio  end  loss  correction  factor  is  given  by 


(9) 


n 


* 


[-*«£>] 


R 


n 


m 


*>(  y~)  I  (1  -  W  ) 


-1 


(10) 


Here  Y\  is  the  normalized  support  temperature,  ^  is  the  measured 
recovery  ratio,  st  and  st*  represent  normalized  resistances,  and  to  is 
a  parameter  which  depends  primarily  on  the  wire  aspect  ratio  and  Nusselt 
number;  CO  represents  the  zero  overheat  value  of  CO  .  The  quantities 
[  t/(t*  -6  )]  and  [  UJ/(  1  -Gd)]  are  given  in  Figures  13  and  14  as 
functions  of  the  two  parameters 


Z?  =  (  £-/d)  f  (k  /k  )  Nu  (1  *  6  a)  ~| 

?  L°W  m(l+st)J 


and 


v  = 


it  S  (£/<o) 

(1  -  O)  ) 
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(11) 


(12) 


It  i8  shown  in  Appendix  A  that  the  correction  factor  iP  ^  does  not  depend 


strongly  on  overheat;  for  f'- 
tractable  form 


0  ,  Eqs.  (9)  -  (12)  assume  the  more 


3*  . 
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( 


*  6 


)  (l^s) 


■  1  +  <  1  *  n 
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s  4  ,  Ca>  , 
}  (  r^io  } 


m 


?  =  (£/d)  f(k0/kw)  Num  (1  +  6  s)  ] 


v  =  -it 


(9a) 

(10a) 

(lla) 

(12a) 


where 


<-.(«»-  >Zm)  <<  i 


t  =  (  n  -  Y)  )  <  <  i 

s  '  1  s  lm' 

and  s  is  a  constant  of  order  one.  As  ^  decreases,  both  and 

depart  from  unity;  the  support  temperature  parameter  v  plays  a  minor 
role  in  determining  the  end  loss  corrections,  although  it  must  be  considered 
if  the  corrections  are  large. 

These  corrections  were  applied  to  the  present  data  and  to  the 


recovery  factor  data  of  Laufer  and  McClellan  at  the  lower  Reynolds 


numbers  of  their  investigation.  As  discussed  in  the  following  section, 
this  correction  was  primarily  responsible  for  unifying  the  results  of  the 
several  investigations  of  recovery  temperature  in  supersonic  flow. 
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For  the  present  data,  the  measured  value  =  0.  903  was  used 

s 

with  the  measured  values  of  /?  and  Nu  to  compute  the  end  loss 

*■  m  m 

correction.  A  slightly  more  complicated  procedure  was  required  in 
interpreting  the  results  of  Laufer  and  McClellan^.  In  that  investigation, 
measurements  of  the  Nusselt  number  and  recovery  factor  were  made  with 
different  wires;  the  Nuq  data  were  obtained  with  wires  of  aspect  ratio 
between  400  and  550  and  the  recovery  temperature  data  were  obtained 
using  aspect  ratios  of  the  order  of  225.  The  Nusselt  numbers  reported 
in  Reference  1  include  a  small  (less  than  5  per  cent)  end  loss  correction; 
for  such  small  corrections  the  formulas  used  are  essentially  identical 
to  those  derived  in  Appendix  A.  The  data  given  in  Reference  1,  therefore, 


represent  Nuq  in  the  present  nomenclature  (the  actual  parameters  used  by 

Laufer  and  McClellan^  were  Nu,  and  Re,;  conversion  to  Nu  and  Re  is 

2  2  o  o 

obvious).  The  reported  recovery  ratio,  however,  represents  anc* 

requires  a  correction  to  account  for  support  conduction. 

The  wire  supports  used  by  Laufer  and  McClellan^  were  equivalent 


to  a  two-dimensional  wedge  of  7  included  angle.  Using  the  observation 


(Appendix  B)  that  the  support  temperature  is  equal  to  the  laminar  boundary 

layer  adiabatic  wall  temperature,  Y\  could  be  calculated  explicitly. 

Since  Nuffi  was  not  known,  an  iterative  computation  process  was  adopted 

which  used  as  input  the  known  quantities  and  the  wire  Reynolds 

number  Re  .  With  Re  known,  Nu  was  taken  from  the  heat  loss 
o  o  o 

correlation  and  the  value  of  Nu^  was  computed.  The  procedure  for 
obtaining  Nu^  and  with  Nuq  and  known  is  analogous  to  the 

determination  of  Nu  and  M.  *with  Nu  and  it  given. 

The  present  free  stream  data  were  subject  to  end  loss  corrections 
in  the  ranges  .703  =  =  .  833  and  1.208  =  =  i-633. 
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For  the  low  Reynolds  number  wake  tests,  values  .  59 


<  5  . 


N 


68 


and  1.70  ;  &  s  2.  06  were  used.  In  view  of  the  appreciable  end 

K. 

losses  encountered,  it  should  be  emphasized  that  the  important  thing  is 
not  the  magnitude  of  the  corrections  per  se,  but  the  accuracy  with  which 
these  corrections  are  known.  The  largest  uncertainty  lies  in  the  actual 
end  conditions  of  the  wire,  as  characterized  by  ;  for  this  reason, 

Y\  was  measured  directly  in  the  present  investigation  (Appendix  B). 

This  parameter  is  sensitive  not  only  to  the  shape  and  material  properties 
of  the  support  but  also  to  the  flow  field  around  the  tips.  In  transonic  or 
non-uniform  flow  fields,  the  determination  of  an  effective  value  of  H 

s 

2 

is  a  formidable  problem,  as  Sherman  has  pointed  out. 

No  correction  was  made  to  the  data  for  radiation  effects  since 
calculation  showed  that  the  error  in  both  Nu  and  Y\  .  was  less  than  one 
tenth  of  one  per  cent.  The  correction  is  roughly  proportional  to  (Tw^/Nu), 
and  may  become  appreciable  for  high  overheats  or  extremely  low  Reynolds 
numbers. 


The  raw  data  from  the  wake  traverse  consisted  of  a  series  of 
continuous  curves  of  wire  voltage  vs.  voltage  output  from  the  traverse 
mechanism.  The  voltage  scale  was  established  by  two  calibration 
measurements  with  the  K-2  potentiometer  at  known  values  of  (y/D). 
Using  this  calibration  and  the  measured  value  of  wire  current,  the  wire 
resistance  was  computed  for  several  representative  stations  within  the 
wake.  Figure  11  shows  the  measured  wire  resistance  as  a  function  of 
wire  current. 

Two  thermodynamic  measurements  were  obtained  with  the  hot¬ 
wire;  these  were  fdR/d(i^)l  .  and  the  zero  current  intercept  R  , 
both  of  which  were  taken  from  Figure  11.  Pitot  pressure  and  static 
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pressure  across  the  wake  were  also  measured  separately.  Thus,  the 
measured  thermodynamic  quantities  were  more  than  sufficient  to 
specify  the  flow  completely,  and  this  redundancy  permitted  an  estimate 
of  the  accuracy  of  the  combined  measurements.  A  description  of  the 
computation  procedure  is  given  below. 

2 

From  the  four  quantities  fdR/d(i  ^  (i.  e.  ,  Nu^), 

Raw  e*  '  Taw*  '  Pt2  '  and  p  ,  four  groups  of  three  independent  quantities 
may  be  formed,  and  any  one  of  these  groups  would  provide  sufficient  data 
to  determine  the  state  of  the  flow  uniquely.  Any  two  groups  would  provide 
an  independent,  redundant  calculation;  in  certain  cases,  the  accuracy 
and  convergence  of  the  iterative  solution  may  depend  strongly  on  the 
choice  of  parameters.  *  For  purposes  of  illustration,  the  following  groups 
of  variables  were  chosen**: 


Group  I:  [  dR/di2J  .=Q  ,  Rawm  ,  p 

Group  II:  p^  ,  Rawm  ,  p 

The  computation  procedure  is  similar  in  principle  to  that  which  would 

be  used  for  any  two  groups  of  three  variables. 

2 

To  interpret  the  local  values  £dR./di  J  in  the  wake  an 
empirical  relation  between  Nuq  and  ReQ  was  required  for  the  range 
.  34  <  ReQ  <  .  66  ,  only  part  of  which  was  covered  by  the  present  free- 
stream  measurements.  For  the  lower  Reynolds  numbers,  the  calibration 
curve  given  in  Figure  8  was  extrapolated  to  ReQ  =  ,  315  using  a  smooth 


*  Several  examples  are  cited  by  Morkovin  to  illustrate  this 

point. 

**  Quantities  without  superscripts  are  functions  of  distance  from 
the  wake  centerline,  while  (  )'  signifies  the  measured  value  at  the  edge 
of  the  wake  mixing  region. 
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curve  which  asymptotically  approached  the  free  molecule  heat  transfer 
law  [Eq.  (4)]  ,  with  a  =  1.0.  The  justification  for  this  choice  of  a  is 
given  in  the  following  section.  Interpretation  of  the  measured  recovery 
temperature  was  made  using  values  of  Yl  *  (Kn)  given  by  the  solid  curve 
of  Figure  10. 

The  static  and  pitot  pressure  measurements  were  obtained  by  John 
F.  McCarthy,  Jr.  of  this  laboratory.  Free  stream  conditions  for  the 
hot-wire  and  pressure  traverses  were 


Hot-wire 

Pressure 


M 

x/l 

00 

5.  74 

9.  0 

5.  74 

9.  0 

Re 


00,  D 


1.  75  x  104 

2.  60  x  104 


Hot-wire  fluctuation  measurements  showed  that  the  wake  was  laminar  in 
both  cases,  and  the  assumption  was  made  that  the  local  pressure  profiles 
were  similar  for  the  two  tests.  The  Reynolds  number  effect  is  discussed 
more  fully  in  Section  V. 

Data  reduction  was  most  conveniently  carried  out  by  assuming 
distributions  of  the  unknown  flow  quantities  across  the  wake,  computing 
the  end  loss  corrections  based  on  these  assumed  values,  and  from  the 
measured  quantities  and  the  end  loss  corrections  obtaining  improved 
guesses  for  the  unknown  quantities.  For  Group  I,  this  iteration  was 
performed  using  M  and  Tq  as  the  unknowns;  the  calculation  for  Group  II 
used  Tq  as  the  iterative  quantity.  In  both  cases,  only  three  iterations 
were  required  to  obtain  convergence. 
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IV.  RESULTS  AND  DISCUSSION 

IV.  1.  Heat  Transfer  Measurements 

In  attempting  to  compile  an  empirical  correlation  of  hot-wire  heat 
transfer  data  over  an  extended  range  of  variables,  the  investigator  is 
faced  with  two  problems:  the  number  of  experimental  papers  on  the 
subject  is  quite  large,  running  well  over  fifty,  and  the  results  of  investi¬ 
gations  conducted  under  seemingly  identical  conditions  differ  by  as  much 

as  40  per  cent.  The  effects  of  wire  overheat  and  stagnation  temperature 
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(particularly  in  the  transonic  regime)  have  yet  to  be  fully  resolved. 

Figure  7  has  been  prepared  as  a  guide  to  the  experimental  variation 

of  Nusselt  number  as  a  function  of  Reynolds  number  and  Mach  number  in 

the  range  0.  2  =  ReQ  =  100.  In  preparing  the  figure  particular  weight  was 

given  to  the  correlation  of  McAdams^  at  M  =  0  and  the  results  of  Laufer 

and  McClellan*  and  the  present  investigation  at  Mach  numbers  greater 

than  two.  Within  the  Reynolds  number  range  1  =  ReQ  =  100,  the  data  of 
21  7 

Baldwin  and  Cybulski  and  Baldwin  are  representative  of  the  information 
available  at  subsonic  velocity.  Spangenberg1  s  data  (Reference  19,  Fig¬ 
ure  17c)  show  a  similar  Mach  number  variation  for  Reynolds  Numbers 

greater  than  one,  and  far  better  agreement  with  References  1  and  20 
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than  the  work  of  Baldwin.  The  cold  wire  data  of  Christiansen  at  low 
Reynolds  numbers  are  also  included.  These  transonic  measurements 
agree  well  with  the  hot-wire  results  for  ReQ  <  1,  and  provide  a  consistent 
set  of  measurements  in  the  near  free  molecule  flow  regime.  At  low 
Reynolds  numbers,  the  correlation  would  be  changed  considerably  if 
incomplete  accommodation  were  present. 

Theoretical  studies  are  of  limited  use  in  preparing  a  heat  transfer 
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correlation.  The  solution  of  Cole  and  Roshko  for  continuum  Oseen 

24 

flow  agrees  well  with  the  experiments  of  Collis  and  Williams  at  low 

25 

Reynolds  numbers.  Levy  has  considered  the  Oseen  flow  problem  with 
temperature-jump  boundary  conditions,  and  the  results  agree  qualitatively 
with  the  experiments  of  Cybulski  and  Baldwin  .  However,  the  velocities 
and  Reynolds  numbers  covered  by  Oseen  theory  are  quite  small. 

Free  molecule  analysis  has  provided  what  appears  from  Eq.  (4) 
to  be  a  straightforward  calibration  for  all  flows  where  the  Knudsen 
number  is  sufficiently  high.  The  fact  that  the  accommodation  coefficient 
a  is  an  arbitrary  parameter  has  been  a  major  weakness  in  all  attempts  to 
apply  this  theory  to  a  given  set  of  data.  *  It  is  not  uncommon  to  find 
investigations  with  an  experimental  scatter  greater  than  30  per  cent  in  which 
the  authors  infer  a  specific  value  of  accommodation  coefficient  by  com¬ 
paring  the  asymptotic  behavior  of  their  Nuq  -  ReQ  curve  to  Eq.  (4),  and 
ascribing  the  difference  to  incomplete  accommodation. 

The  heat  transfer  data  obtained  in  this  investigation  are  presented 
in  Figure  8.  These  data  agree  well  with  the  measurements  of  Laufer 
and  McClellan^  at  Reynolds  numbers  ReQ  greater  than  three.  One  may 
infer  from  this  agreement  and  the  Mach  number  independence  principle 
that  the  solid  line  shown  in  Figure  8  represents  the  Nuq  -  ReQ  relation 
for  Mach  numbers  greater  than  about  two. 

For  comparative  purposes,  the  results  of  free  molecule  theory 
[Eq.  {4)j  and  the  extrapolated  slope  of  the  high  Reynolds  number  data 

are  shown  in  Figure  8.  At  high  Reynolds  numbers,  the  Nusselt  number 

i. 

becomes  proportional  to  (Reo)2,while  Nuq  ^  ReQ  for  ReQ  <  <  1.  The 

*  An  excellent  appraisal  of  the  problems  of  momentum  and 
energy  accommodation  is  given  by  Hurlbut^. 
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present  measurements  asymptotically  approach  the  free  molecule  result 

with  an  accommodation  coefficient  a  near  unity.  For  the  heat  transfer 

23  6 

tests  in  air  and  argon  by  Christiansen  ,  Weltman  and  Kuhns  , 

3 

and  Wong  ,  an  accommodation  coefficient  of  unity  represents  the  data 

within  the  experimental  error.  Experimentally,  it  has  been  found  that 

surface  conditions  are  of  primary  importance  in  determining  the  degree 
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of  energy  accommodation  ;  the  presence  of  adsorbed  gases  on  the  body 

surface  increases  the  accommodation  coefficient. 

All  data  shown  in  Figure  8  are  based  on  the  measured  Nusselt 

number  at  zero  overheat;  however,  the  effects  of  overheat  are  discussed 
in  detail  in  Appendix  C.  In  free  molecule  flow,  the  surface  tem¬ 


perature  enters  the  heat  transfer  relation  only  through  its  effect  on 

the  accommodation  coefficient.  If  the  accommodation  coefficient  were 

unaffected  by  surface  temperature,  the  cold  wire  results  of  Christiansen 

and  the  present  results  should  agree  closely  [except  for  the  small  Mach 

number  variation  given  by  f(Sj)/sj  in  Eq.  (4)[j  .  Figure  7  indicates 

that  the  hot  wire  and  cold  wire  data  agree  well  up  to  a  Reynolds  number 

Re  of  about  one. 
o 

For  high  Reynolds  number  flows,  a  large  part  of  the  heat  transfer 

to  the  body  occurs  in  the  vicinity  of  the  forward  stagnation  point.  For  a 

perfect  gas,  the  heat  transfer  at  the  stagnation  point  of  a  cylinder  in 

27 

supersonic  flow  is  closely  represented  by 
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Equation  (13)  indicates  that  the  effects  of  overheat  at  high  Reynolds 
numbers  should  be  small  over  the  range  of  conditions  encountered  in  the 


laboratory. 
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Although  no  complete  solutions  are  available  for  a  cylinder  in 
high  speed  flow  because  of  the  lack  of  an  adequate  theory  covering  the 
region  of  separation,  many  analyses  are  available  for  heat  transfer  near 
the  stagnation  point  of  a  blunt  body.  It  has  been  pointed  out  that  both 
the  hot-wire  data  and  the  stagnation  point  boundary  layer  solutions 
exhibit  heat  transfer  coefficients  proportional  to  /Re  at  high  Reynolds 
numbers,  and  the  behavior  of  the  stagnation  point  solutions  with 
decreasing  Reynolds  number  should  agree  well  with  the  hot  wire 
measurements. 

3  2 

Two  recent  investigations  by  Probstein  and  Kemp  and  Van  Dyke 
have  used  the  Navier -Stokes  equations  to  describe  the  flow  field  near  the 
stagnation  point  of  a  blunt  body  in  hypersonic  flow.  These  analyses  show 
that  conventional  high  Reynolds  number  boundary  layer  theory  must  be 
modified  as  the  body  Reynolds  number  is  decreased.  First  order  corrections 
may  be  present  due  to  external  vorticity,  boundary  layer  displacement, 
surface  velocity  and  temperature  jump,  and  body  curvature.  For  a 
cylinder,  the  vorticity  along  the  zero  streamline  is  zero  in  the  inviscid 
flow,  so  the  vorticity  effect  on  heat  transfer  is  of  second  order.  However, 
the  boundary  layer  displacement  and  curvature  effects  both  produce  a 
first  order  decrease  in  surface  heat  transfer,  compared  with  the  con¬ 
ventional  high  Reynolds  number  value.  Slip  effects  would  also  decrease 
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the  heat  transfer  rate.  Probstein  and  Kemp  have  proposed  a  merged 
layer  analysis  based  on  continuum  relations  for  the  range  30  -S  Re  ^/100. 


The  results  agree  qualitatively  with  the  present  measurements. 
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The  slope  of  the  Nuq  -  Re^  relation  is  a  quantity  of  considerable 
interest.  The  exponent  n  in  the  relation 

Nu  ~  Re  n  (14) 

o  o  ' 

correlates  the  results  of  several  independent  investigations  better  than 

the  absolute  magnitudes  of  the  observed  quantities.  In  Figure  9,  the 

exponent  n  is  shown  as  a  function  of  Reynolds  number  and  Mach  number 

for  several  investigations.  The  present  measurements  and  those  of 

Eaufer  and  McClellan^  were  used  for  M  >  2,  while  selected  results  of 

.  .  22  7  21 

Christiansen  ,  Cybulski  ,  and  Baldwin  are  shown  for  M  <  1.  At  high 

Reynolds  numbers,  the  exponent  n  is  in  the  limit  of  Kn  — 0,  the 

free  molecule  relation,  Eq.  (4),  shows  that  n  =  1.  A  third  relation, 

23 

derived  from  the  Oseen  solution  of  Cole  and  Roshko  ,  is  applicable 
when  M/Re  <  <  1  and  ReQ  <  <  1;  in  this  instance,  the  exponent  is  of  the 
form 


n  = 


log  ( 


(15) 


Re  Pr 
o  o 
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Collis  and  Williams'  experiments  agree  very  well  with  Eq.  (15)  for 


0.  01  <  Re  <  0.  3. 
o 

At  high  Mach  numbers,  the  exponent  n  is  seen  to  vary  monotonically 

between  the  continuum  and  free  molecule  values  of  \  and  1.  The  measured 

slope  begins  to  depart  from  the  continuum  value  at  a  Reynolds  number 

ReQ  of  about  200  for  an  insulated  wire.  A  comparison  between 

22 

Christiansen's  data  and  the  present  measurements  indicate  that  the 
departure  from  the  continuum  heat  transfer  law  occurs  at  a  slightly 
lower  value  of  ReQ  if  the  wall  is  highly  cooled;  this  may  be  explained 
by  the  fact  that  a  cool  wall  decreases  the  boundary  layer  displacement 
effect  at  a  given  Reynolds  number.  For  low  Reynolds  numbers,  the 
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approach  to  the  free  molecule  heat  transfer  law  does  not  appear  to  vary 

with  wire  overheat,  although  the  data  are  not  sufficient  to  completely 

resolve  this  question. 
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Willis  has  treated  several  problems  in  the  near  free  molecule 
regime  using  an  iterative  solution  to  the  integral  form  of  the  Boltzmann 
equation.  This  approach  is  equivalent  to  the  so-called  first  collision 
theory  if  only  the  first  iterate  is  considered.  Willis'  results  for  a 
two-dimensional  strip  normal  to  a  rarefied  flow  should  be  identical  in 
form  to  the  solution  for  a  cylinder  normal  to  the  flow.  By  writing  a  heat 
balance  for  the  body,  the  resulting  equations  may  be  solved  for  the 
departure  of  the  Nusselt  number  and  recovery  temperature  from  the 
free  molecule  values  as  a  function  of  the  free  stream  Knudsen  number. 

The  results  may  be  expressed  as  (a  =  1) 

Nu 

°  =  1  +/aX  (16) 

ee  Molecule 


(17) 

(18) 

(19) 


Here  is  the  incident  molecular  flux,  is  the  incident  energy  flux, 
and  (  indicates  the  value  for  LL  — >  0  . 


(20) 
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Willis  has  evaluated  the  parameters  and  %  f°r  a  tw°* 

dimensional  strip,  and  finds  that  they  depend  on  the  molecular  speed 

ratio,  the  body  temperature,  the  molecular  model,  and  the  number  of 

degrees  of  freedom  of  the  gas.  and  X.  are  weak  functions  of  the 

body  temperature  if  the  overheat  is  small;  the  results  are  more  sensitive 

to  body  temperature  at  higher  values  of  the  molecular  speed  ratio  Sj. 

For  Sj  — ^  oo  ,  Willis  has  computed  ^  and  ~X.  as  a  function  of  body 

temperature  for  both  the  hard  sphere  and  modified  Krook  molecular 

models.  For  an  insulated  body  (Tw  =  Taw)  ,  ^  =  .  42  and  X  =  •  11  for 

the  modified  Krook  model  while  ^  =  .  78  and  X  =  .  48  if  the  hard 

sphere  model  is  used.  From  Equation  (16),  it  may  be  seen  that  the 

departure  of  the  Nusselt  number  from  the  free  molecule  result  is  a 

strong  function  of  the  molecular  model  employed.  Eqs.  (16)  and  (17) 

agree  qualitatively  with  the  present  experiments  if  ^  >  X  >  0  • 

Several  interesting  conclusions  may  be  drawn  from  the  above 

analysis.  First,  the  small  parameter  in  near  free  molecule  flow  for  a 

two-dimensional  body  is  (l/Kn)  log  (l/Kn)  rather  than  the  quantity 

(l/Kn)  which  pertains  to  axi symmetric  bodies.  Thus,  the  range  of 

validity  of  the  two-dimensional  analysis  is  restricted  to  higher  Knudsen 

numbers  than  the  axi- symmetric  counterpart.  Second,  for  viscosity 

1 

proportional  to  (T)2  ,  Kn  ^  (l/ReQ)  at  high  Mach  number  and  the 
analysis  predicts  that  the  Nuq  -  ReQ  relation  at  high  Knudsen  numbers 
should,  as  expected,  become  independent  of  Mach  number. 
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IV.  2.  Recovery  Temperature  Measurements 

Measurements  of  the  recovery  temperature  of  a  cylinder 
placed  normal  to  the  flow  have  previously  been  made  at  high  Reynolds 
numbers  up  to  a  Mach  number  of  4.  5.  The  excellent  work  of  Laufer  and 
McClellan  has  shown  that  the  recovery  temperature  is  independent  of 
Mach  number  above  approximately  2.  A  justification  of  this  effect  has 
already  been  given. 

In  the  transitional  regime,  the  measurements  reported  in  References 
1  and  2  indicated  that  at  Knudsen  numbers  of  about  0.  1  the  recovery  tem¬ 
perature  begins  to  rise  above  the  high  Reynolds  number  value,  while  at 
Knudsen  numbers  greater  than  10,  the  full  temperature  rise  predicted  by  Eq. 
(5)  is  attained.  However,  at  Knudsen  numbers  of  about  1  and  at  roughly  the 
same  Mach  number,  the  two  investigations  reported  significant  differences 
in  the  recovery  factor. 

The  present  investigation  was  undertaken  in  an  attempt  to  resolve 
this  apparent  discrepancy  and  at  the  same  time  provide  additional 
recovery  temperature  data  for  use  in  steady- state  hot  wire  work.  The 
results  are  shown  in  Figure  10,  using  the  normalized  recovery  ratio  ^ 
and  the  free- stream  Knudsen  number  as  parameters.  Fortunately,  the 
recovery  temperature  for  a  cylinder  is  independent  of  the  accommodation 
coefficient  (provided  that  the  wire  is  infinitely  long.')  and  rapidly  approaches 
an  asymptote  as  the  Mach  number  increases.  Therefore,  since  Y\  ^  =  0 
for  H  and  Y\^  =  1  for  =  Y\^  ,  all  investigations 

should  agree  at  least  in  these  two  limits. 

The  shaded  area  representing  the  data  of  Laufer  and  McClellan 
defines  the  limits  of  experimental  scatter  of  nearly  100  separate 
measurements  made  over  a  large  range  of  Mach  and  Reynolds  numbers. 
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Considering  the  range  of  parameters  studied,  the  data  are  extremely 
consistent.  However,  the  aspect  ratios  of  the  wires  used  in  their 
recovery  temperature  study  were  on  the  order  of  200,  and  at  the  higher 
Knudsen  numbers  these  wires  were  subject  to  significant  end  losses. 

Using  the  relations  derived  in  Appendix  A  several  representative  data 
points  were  corrected  for  support  conduction  and  plotted  in  Figure  10. 

By  combining  the  corrected  data  of  Reference  1,  the  results  of 
2 

Sherman  ,  and  the  present  measurements,  a  continuous  transition  of  the 
normalized  recovery  factor  is  obtained  between  continuum  and  free 
molecule  flow.  From  a  consideration  of  the  Mach  number  independence 
principle,  this  curve  should  be  valid  for  all  Mach  numbers  above 
approximately  two.  Considering  the  sensitivity  of  the  normalized  recovery 
ratio  to  experimental  error,  the  correlation  appears  to  be  very  good. 

At  lower  Mach  numbers,  the  ft*  (Kn)  correlation  should  show  a 
dependence  on  Mach  number. 

The  behavior  of  the  *(Kn)  relation  at  small  Knudsen  numbers 
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agrees  qualitatively  with  the  results  of  Probstein  and  Kemp  .  In  the 
incipient  merged  layer  analysis,  they  find  that  the  gas  temperature  at 
the  stagnation  point  of  a  blunt  body  becomes  greater  than  the  free  stream 
stagnation  temperature  in  the  low  Reynolds  number  continuum  range. 

For  large  Knudsen  numbers,  Eq.  (17)  predicts  that  the  recovery  tem¬ 
perature  will  fall  below  the  free  molecule  value  if  >  X  l  the 

experimental  data  shown  in  Figure  10  indicate  that  for  a  circular 
cylinder  and  a  diatomic  gas,  ^  and  X  are  nearly  equal. 
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V.  APPLICATION  OF  THE  STEADY-STATE  HOT-WIRE 
TO  WAKE  MEASUREMENTS 


Several  classic  experiments  have  been  made  using  the  hot-wire 
as  a  steady- state  measuring  instrument.  These  have  been  well  reported 
in  the  cited  references,  and  it  is  pertinent  only  to  emphasize  that  wire 
calibration  is  the  primary  problem  in  all  quantitative  measurements. 

The  number  of  investigations  utilizing  the  hot-wire  in  a  redundant 
measurement  scheme  is  quite  small,  and  consequently  it  was  decided  to 
perform  an  illustrative  experiment  involving  four  measured  quantities, 
two  of  which  were  obtained  from  a  hot-wire  operated  in  the  transitional 
regime.  The  wake  of  a  cylinder  transverse  to  a  hypersonic  flow  was 
chosen  both  because  its  flow  field  was  amenable  to  hot-wire  and  pressure 
measurements  and  because  of  current  interest  in  the  thermal  wakes  of 
hypersonic  blunt  bodies. 

The  wire  used  in  this  experiment  had  an  aspect  ratio  of  311  and 

2 

a  =0.  93.  The  measured  variation  of  wire  resistance  with  i  is  shown 
r 

in  Figure  11.  One  important  conclusion  which  can  be  drawn  is  that  the 

instrument  sensitivity  was  quite  sufficient  to  pick  out  small  differences 
2 

in  the  R  vs.  i  relation.  This  property  is  absolutely  essential  for 
quantitative  interpretation. 

The  following  groups  of  variables  were  chosen  to  represent  the 


data.-* 

Group  I: 
Group  II: 


dR 

-d(i2)  - 

Pt2  ’  P  ' 


,  R 

awm 

L  =  o 

R 

awm 


P 


*  See  Section  III  for  a  discussion  of  the  data  reduction  procedure. 
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The  distributions  of  flow  quantities  across  the  wake  were  computed 
using  the  measurements  of  Group  I  and  Group  II  in  independent  calcula¬ 
tions.  Thus,  the  agreement  between  the  results  of  Group  I  and  Group  II 
provides  a  measure  of  the  accuracy  of  the  combined  measurements. 

Figure  12a  shows  the  results  of  the  stagnation  temperature 
calculation.  T  1  is  the  measured  stagnation  temperature  at  the  outer 
edge  of  the  mixing  layer.  Results  for  the  several  representative  stations 
are  shown  first  as  "uncorrected"  data.  By  "uncorrected"  is  meant  the 
temperature  one  would  calculate  using  only  the  zero  current  intercept 
of  Figure  11  and  the  wire  resistivity  coefficient.  The  wire  was  subject 
to  end  loss  effects,  and  by  using  Figures  8  and  10  and  the  derived  end 
loss  corrections,  an  iteration  was  performed  to  find  the  true  value  of 
Tq.  Only  three  iterations  were  required  to  obtain  convergence. 

Qualitatively,  the  behavior  of  (Tq/T  ')  is  similar  to  that  obtained 
in  a  high  speed  boundary  layer  *  with  Prandtl  number  less  than  one. 
Because  of  laminar  shear  and  heat  conduction  processes,  the  centerline 
Mach  number  for  this  particular  streamwise  station  was  about  3/4  of  M’ 
and  the  centerline  stagnation  temperature  was  about  0.  96  of  T  The 
two  values  of  T  1  agreed  within  0.  14  per  cent,  and  were  two  per  cent 
lower  than  the  free  stream  stagnation  temperature.  This  result  agrees 
well  with  the  postulate  of  isoenergetic  flow  everywhere  outside  the 
viscous  mixing  region. 

The  local  mass  flow  in  the  wake  is  shown  in  Figure  12b.  For 

purposes  of  comparison,  the  "uncorrected"  hot-wire  data  is  also  shown; 

/  2  -1 

the  measured  heat  loss  of  the  wire  R&w  (dR/di  )  is  directly  proportional 
to  the  local  mass  flow  (pu)  if  end  losses  are  negligible.  Since  the  end 
loss  correction  increases  with  decreasing  (pu),  the  uncorrected  mass  flow 
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is  higher  than  the  final  corrected  values  based  on  the  measurements  of 
Group  I. 

The  agreement  between  the  two  values  of  (pu/p'u*)  computed 

using  Group  I  and  Group  II  is  good.  The  differences  between  the  two 

curves  are  primarily  due  to  the  different  cylinder  Reynolds  numbers 

for  the  hot-wire  and  pressure  tests  (see  page  20).  Subsequent  pitot  pressure 

data  have  indicated  that  the  higher  value  of  (pu/p'u1)  at  the  wake  centerline 

obtained  with  Group  II  is  caused  by  the  higher  values  of  Re  for  the  pres- 

00,1 

sure  measurements.  At  the  edge  of  the  wake,  (p'u1)  for  Group  I  was  0.  85  of 

(p'u1)  for  Group  II.  Taking  this  Reynolds  number  effect  into  account,  the 

absolute  values  of  (  pu)  from  Group  I  and  Group  II  were  found  to  agree 

within  7  per  cent  throughout  the  mixing  region.  The  small  differences 

between  the  two  values  of  (T0/To*)  shown  in  Figure  12a  are  also  attributable 

to  the  difference  in  Re  ^ 

oo,  D 

The  velocity  defect  in  the  wake  mixing  region  is  shown  in  Figure 

12c.  The  agreement  between  the  hot-wire  and  pressure  measurements 

is  seen  to  be  excellent,  and  verifies  the  fact  that  the  differences  shown 

in  Figures  12a  and  12b  are  due  to  differences  in  Re  ^  ;  for  (x/D)  =  9 

oo,  -LI 

and  these  Reynolds  numbers,  the  wake  is  laminar  and  the  velocity  defect 
[(u  -  u<£  )  /  (u1  -  u^  )J  should  show  little  dependence  on  Re^  -p  • 

In  summary,  this  wake  survey  has  shown  that  the  steady- state 
hot-wire  may  be  used  as  an  accurate  quantitative  instrument  when 
operated  at  small  wire  Reynolds  numbers,  provided  sufficiently  accurate 
calibration  relations  are  available.  The  hot-wire  is  capable  of  providing 
two  thermodynamic  measurements;  only  one  additional  measurement 
(such  as  total  or  static  pressure)  is  required  to  determine  the  flow  field 
uniquely. 
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APPENDIX  A 

END  LOSS  CORRECTIONS  FOR  THE  TRANSITIONAL  REGIME 

In  interpreting  hot-wire  measurements  in  low  density  flows,  it 
is  imperative  that  the  heat  loss  to  the  wire  supports  be  accurately  taken 
into  account.  For  high  Reynolds  number  flow  where  the  convective  heat 
transfer  is  large  and  end  loss  corrections  are  small,  the  derivation  given 
by  Kovasznay^  is  adequate.  However,  considerable  care  must  be  exer¬ 
cised  when  these  corrections  amount  to  20  per  cent  or  more  of  the 
measured  quantity. 

Formulation  of  the  Problem 

Neglecting  radiation,  we  may  write  the  following  heat  balance 
for  a  small  wire  element  of  length  dz: 

»d^h(T-Tawni  =  ){d/dz)  (kw  «  )  +  i2  R  (A-I) 

where  the  left  side  of  the  equation  represents  the  convective  heat  loss, 

the  first  term  on  the  right  hand  side  is  the  net  heat  added  by  conduction, 

2 

and  (i  R)  is  the  electrical  dissipation.  Solution  of  this  equation  involves 
the  following  assumptions: 

(a)  h,  the  convective  heat  transfer  coefficient,  is  a  constant  and 
is  equal  to  the  value  which  would  be  measured  by  an  infinite  wire  at  the 
same  current.  This  implies  that  the  external  flow  is  uniform  along  the 
wire . 

(b)  The  wire  thermal  conductivity,  kw  ,  is  constant. 


3? 


(c)  Wire  resistance  varies  linearly  with  temperature: 
R 


1  +  a  (T  -  T  ) 
r  '  r 


=  1  +  st 


'’awm  1  +  a  (T  -  T  ) 
r  '  awm  r' 


1  +  ar  (Tawm  "  Tr^ 


(A- 2) 


{T  -  T  )/T 
'  awm''  o 


(d)  Each  cross-section  of  the  wire  is  at  uniform  temperature 

£  implicit  in  Eq.  (A-l)  ]  . 

(e)  Both  wire  supports  are  maintained  at  a  temperature 
Tg  =  Y\  Tq  .  *  The  boundary  conditions  are  thus 

z=  -  j  ,  t=  tw=  (Tg  -  Tawm)/TQ 


(A-3) 


With  the  above  boundary  conditions,  Eq.  (A-l)  may  be  integrated  to  give 

k 


t  = 


P  s  w  COS  Lily 


+  t  cosh  2^  (z/j!  ) 
s  cosK  v 


where  2 

2  t  £  \  ^  *  Rawm  ^  S 

v  --  (|)  <ir>Nuo - yr—z - 

. .  ^  w  ird  k  T 

e  =  w  0 

The  average  wire  temperature  at  any  given  current  is 


(A-4) 


JL 


1/z 

k 

J  tdz  ’  -Z- 1 

<f>2 

{  ^  )  nuq  (i  +  e  s) 

2]  |  tanh  V  ~| 

- »  J  L1  -  >  J 

-A/2  v  3 

w 

4-  f 

tanh 

+  s 

(A-5) 


*  This  restriction  is  not  essential  although  it  simplifies  the 
mathematics  considerably.  The  effect  of  non-uniform  flow  near  the 
supports  may  be  partially  accounted  for  by  an  appropriate  choice  of 
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For  an  infinitely  long  wire,  {  £,  /d)  — ►  oo  ,  t  — ^  t+  ,  [(tanh  V/ u  5 — ►  0 
and 

k 

Um  r  2  t„  =  -V  [  (  f  >2  <  )Nuo(l  +  fs)-  V2]  . 

(£,/d)-f  oo  P  s  u  w 

Eq.  (A-5)  may  now  be  rewritten  in  the  final  form 


7-  _  t  ,  ,  tanh  V  4  ,  4  tanh  >* 

1  -  \  1  “  “  /  +  lt 


(A- 6) 


18  ,, 


JS  '  s  j/ 

This  relation  is  identical  to  that  derived  by  KovasznayiU  if  we  let  the 

support  temperature  be  equal  to  the  wire  recovery  temperature  (t  =  0) 

and  £  =  0.  It  may  easily  be  seen  that  the  wire  temperature  t  will  always 

lie  between  the  temperature  t^  of  an  infinite  wire  and  the  support 

temperature  t  .  For  large  end  losses  (  V  small),  the  mean  wire 
s 

temperature  approaches  the  temperature  of  the  supports. 


Nusselt  Number  Correction 


In  order  to  correct  Nu  to  the  value  Nu  for  an  infinite  wire,  we 

m  o 

must  find  the  ratio  of  these  two  quantities  at  a  given  current.  In  the 
present  nomenclature,  this  ratio  is 


f 


Nu 


N 


Nu 


=  ( 


m 


Introducing  the  quantities 


(1  +  st*) 
(1  +  sT  ) 


(A- 7) 


(  v  )i=Q  =  V(1  +  St  )»  f 
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and  defining 


CO  -  (tanh  !>  /  D  )  ,  CO  -  (tanh  V  /  £?  ) 
we  have 

•V 

6  »  -  t  <  )  .  (A- 8) 

s  1  -  CO 

Equations  (A-6)  and  (A- 8)  must  be  solved  simultaneously  for  t^  and 
in  terms  of  the  measured  quantities  of  the  experiment. 

For  low  overheats  ,  we  may  perform  an  expansion  in  the 
parameter  (st"  ).  From  Eqs.  (A-6)  and  (A-8), 


(  ——"V  )  =  (1  -  ^  )  [  1  +  G  (  f  ,  t  ,  tfl  )] 


where 


G 


(1  -CO)[l  -  Co) 


Using  the  relations 


and 


(?/»)  I  st  -  J(sF)2  + 


tanh  (  V  -  J J  ) 


tanh  l>  -  tanh  U 
1  -  tanh  V  tanh 


1 


a  certain  amount  of  manipulation  leads  to  the  equation 


(A-9) 


(A- 10) 


G  =  its 


CO 


(1  -W)(l  -CO) 


X 


1 11  (s  t)  +  ••• +  —  kr^l 

...)] . 


I- 1  _  (1  -  tanh  l>  tanh  V  )  f  ^  V  .  Z 

L  “  £  1  “T_  » 


x 

(A- 11) 


At  this  point,  we  are  able  to  make  some  definite  arguments  about  the 
magnitudes  of  the  terms  appearing  in  Eq.  (A-9).  The  quantity  s  will 
in  general  be  of  order  1  and  if  the  end  loss  is  to  be  small,  then  Co 


must  also  be  small  and 
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<  <  1. 


(  J  -  <o)(l  -CO) 


Since  |tgj  =  J  (  ^  s  “  ^  m)  f  <<  1  ,  G  will  then  be  °  (  ft  m  "  K  s> 
if  the  bracketed  quantities  in  Eq.  (A-ll)  are  0(1).  The  following 
approximations  are  now  made: 

(a)  The  quantity  (  -  1 )/(t>j(  -  £  )  is  small,  justified  by  the 

experimental  fact  that  Nusselt  number  changes  slowly  with  overheat. 

(b)  (1  -  tanh  5?  tanh  V  )  <  <  1.  This  is  strictly  true  only  for 
large  ,  but  making  this  assumption  places  an  upper  bound  on  G  and 
may  be  checked  a  posteriori. 


Defining 


v  =  - 


ft  s  v  (£/co  )  SLz 
s  o  (1  - 


(i  -co) 

(1  -co  ) 


(A- 12) 


we  have  the  form 


v  co 


( 1  -  tanh  tanh  >}]  (A-13) 

Substituting  Eq.  (A-13)  into  Eq.  (A- 9)  and  noting  that  v  <  <  1,  we  obtain 


the  final  relations 


N  Nu. 


■  <  T 


(i  +  St*) 

(  1  +  St) 


(A- 14) 


(1  -60)  [l - —  Z-  U  + 

L  ( 1  -60)  1 


(A- 15) 


.  .  .  +o((st  )  ,  J-f >  (i  ~  tanh ))  tanh  v5  )  ,  v2^}] 


=  ~  I  tg  s  |  (  td/ to  ) 


(1  -co) 


(A- 16) 


>¥< 


(A- 17) 
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=  U/d) 

r  (k  /k  )  Nu 

'  o  w  m 

2 

(1  +  6  s) 

L 

-  (1  +  s  T  )  J 

I2  ■  [<  TT-7  >  Nv]  i=o  [«  W  >  NV 


__  rv 

For  the  special  case  of  zero  overheat,  t  =  0  ,  t^  =  e  ,  V  =  >>  , 
y  =  1  ,  and  Eqs.  (A-6)  and  (A-8)  may  be  combined  to  give  an  explicit 


form  for  the  Nusselt  number  correction  factor: 


<^N>i=o  *  <l  '")[ 


i  -  it  s  ( 


CO 


)  -  4  (t 

CO  '  2  '  S 


S)  ( 


CO 

i  -  co 


(A- 20) 


Recovery  Factor  Correction 

From  Eq.(A-8),  we  may  easily  find  the  recovery  factor  corrections 


to  be 


f  , 

K  * 

r  i  -  <  n  j 

1 

rR  ~ 

^2-  m 

1  (i  -  <2  ) 

J 

- 

K  *  _ 

,+  .a  r,  _  . 

nc- 

Ks 

R 

K 

u  m 

(1  -U»  L 

Yl  rn  - 

nc 

(A-21) 

(A-22) 


Computational  Procedure 


At  this  juncture,  it  seems  advantageous  to  illustrate  the  com¬ 
putational  advantages  of  Eq.  (A- 15).  Neglecting  higher  order  terms, 
if  v  and  V  are  taken  to  be  parameters  we  have  the  relations 


f  (  V,  v) 


t*  - 


V  =  f 


6 
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and  hence 

<  -rr  i  ■  f  <<■  v  i 

V  =  f  <£,  V  )  , 

In  Figures  13  and  14  ,  [  t /(t^  -  £  )]  and  [OJ /( 1  -  <^>)]  are  shown  as 
functions  of  (l/ ),  with  v  as  a  parameter.  The  quantity  S  used  by 
Kovasznay  is  related  to  by 

S  =  lim  (1/^  )  .  (A- 23) 

£  -vO 

The  analysis  of  that  reference  further  restricts  the  solution  to  cases 
where  ^  =  1  and  tg  =  0. 

It  may  be  readily  seen  from  Figure  13  that  the  parameter  v 
becomes  important  only  in  circumstances  where  the  end  loss  correction 
is  sizeable.  To  assess  the  relative  magnitude  of  v  and  its  effect  on  the 
solution  we  take  a  hypothetical  case  where  the  end  loss  may  be  expected 
to  be  large.  Let 

(k  /k  ) 2  =  03125,  Nu  =  0.  1600  ,  t  =  0  ,  £  =  1 

'  o'  w'  '  m  *  ’  j 

t  =  -0.  120  ,  s=f,  (£/d)  =  250 
s 

nc  =  .950  ,  }fm=  1.020  ,  >?f-  >lc  =  0.210  . 

To  start  the  calculation,  assume  6=0.  Then  from  Eq.  (A- 18) 

3T  ^  =  .03125  x  0.4  x  250  =  3.125 
=  . 3200  . 

From  Eq.  (A- 16),  we  calculate 

v  =  -  \  (-0.  120)  .  |  =  .  030 
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and  from  Figures  13  and  14  we  obtain 

<r/t„  -6)<0>  =  0.571  .  [  '('r'Ho )  ]  <0)  =  °'710  • 

Using  Eq.  (A- 15)  with  t  ^  0  and  t^  =  £  =  0  , 

N<0>  =  0.571  . 

The  first  approximation  to  the  correction  factor  for  Y[  ^  is  computed 
from  Eq.  (A-22): 

^R<0)  *  I  V  'Ll101  =  1  +  0.710  [l+  =  2.217 

and 

e<0>  =  («f  -  nc)  «m<  -2-1  -n 

^  m 

=  0.  210  x  0.  333  x  1.  217 
=  .0852 

The  calculation  may  now  be  repeated  and  a  second  approximation  obtained 
using  The  following  results  are  obtained: 


Quantity 

Assumed  Value 

(  >(0) 

<  >“> 

<  ><2) 

6 

0 

0.  0852 

0.  0817 

0.  0817 

fR 

- 

2.  217 

2.  167 

2.  167 

TP,, 

- 

0.  571 

0.  608 

0.  608 

N 


This  hypothetical  example  illustrates  several  general  properties 
of  the  solution.  Even  for  large  corrections,  the  iteration  scheme 
rapidly  converges  to  the  true  value,  regardless  of  the  initial  choice  of 
6  .  Further,  the  magnitudes  of  the  neglected  terms  in  Eq.  (A-15)  may 
be  evaluated  and  suitable  correction  applied  before  the  final  iteration. 
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For  finite  overheats,  ^  will  also  enter  the  computation;  since  $  appears 
only  as  a  multiplier  of  v  and  is  always  close  to  one,  its  contribution  to 
the  computation  should  be  negligible. 

Summary 

A  rapidly  convergent  procedure  for  obtaining  end  loss  corrections 
appropriate  to  transitional  flow  is  given.  Eq.  (A-22)  indicates  that  the 
support  temperature  may  contribute  significantly  to  the  correction  of  the 
normalized  recovery  ratio,  and  unless  the  end  losses  are  small,  it 
must  be  accurately  known.  The  support  temperature  enters  the  Nusselt 
number  correction  only  through  the  parameter  v,  and  unless  the 
parameter  is  less  than  5,  it  has  a  negligible  effect  on  the  corrected 
value  Nu 
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APPENDIX  B 

MEASUREMENT  OF  THE  WIRE  SUPPORT  TEMPERATURE 

The  analysis  given  in  Appendix  A  shows  clearly  that  the  temper¬ 
ature  of  the  hot-wire  supports  is  an  important  parameter  in  the 
interpretation  of  recovery  temperature  data.  This  is  particularly  true 
when  the  normalized  parameter  ^  is  used,  since  an  error  in  the 
determination  of  the  absolute  wire  temperature  is  magnified  by  a  factor 
of  five  or  more  in  the  computation  of  Y\  ^  .  For  this  reason,  a  special 

test  was  conducted  to  determine  the  quantity  Y\  explicitly, 

s 

Two  .  001"  diameter  wires,  one  Pt  -  10  per  cent  Rh  and  one 
Constantan,  were  soft- solder ed  to  the  tip  of  one  of  the  hot-wire  support 
needles  to  form  a  thermocouple  junction.  The  complete  assembly  was 
then  placed  in  the  calibrating  oven  and  a  curve  of  e.  m.  f.  vs.  junction 
temperature  obtained.  The  thermocouple  wires  were  terminated  in  a 
constant- temperature  junction  box  for  both  the  oven  and  tunnel  tests, 
thus  eliminating  spurious  e.  m.  f,  sources  caused  by  thermocouple  wire- 
lead  wire  temperature  variation. 

The  hot-wire  support  and  thermocouple  combination  was  placed 
in  the  tunnel  to  reproduce,  as  closely  as  possible,  the  actual  flow  con¬ 
ditions  encountered  in  the  hot-wire  measurements.  Since  the  needle 
supports  tapered  from  ,  02”  diameter  at  several  calibers  to  about  .  01" 
at  the  tip,  the  .  001"  thermocouple  wire  was  small  compared  to  body 
dimensions  and  should  have  had  a  small  effect  on  the  surrounding  flow 
field.  No  correction  was  made  for  heat  loss  by  conduction  along  the 
thermocouple  wire. 
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Measurements  of  the  support  temperature  were  made  over  the 
complete  range  of  tunnel  conditions  for  several  positions  of  the  thermo¬ 
couple  near  the  needle  stagnation  point.  These  data  are  shown  in 
Figure  15.  The  needle  tip  was  roughly  ogival  in  shape  (as  shown  in  the 
insert)  and  the  angle  9  used  to  define  the  thermocouple  location  refers 
to  the  angle  between  the  surface  normal  and  the  free- stream  direction 
at  that  point.  These  data  indicate  that  the  support  temperature  was  a 
constant  and  equal  to  0.  903  Tq  over  the  complete  operating  range. 

This  result  can  be  corroborated  by  a  simple  dimensional 
argument.  The  rate  at  which  heat  is  supplied  to  the  probe  by  convection 


is  proportional  to 


- 


k  Nu  (AT  ) 
a  o'  g' 


■1 


where  A  T^  is  the  local 


difference  between  the  actual  support  temperature  and  the  adiabatic 
temperature.  From  point  to  point  within  the  solid,  the  heat  transfer 


rate  is  proportional  to 


k  (  AT  ) 


x 


where  A  T  is  the  temperature 


difference  between  two  points  a  distance  x  apart  and  ko  is  the  thermal 
conductivity  of  the  support.  The  ratio  of  the  two  temperature  differences 
is  thus  of  order 


(ATs/  A  Tg)  ~  (ka/ks)  (x/d)  Nuo  . 

The  ratio  (k  /k  )  for  air  and  a  steel  needle  is  of  order  10  ^  ;  hence 

cl  S 

AT  <  <  AT  ,  and  the  needle  is  very  close  to  constant  temperature 

s  S 

throughout. 

Let  us  carry  the  argument  one  step  further.  The  free-stream 

3 

Reynolds  number  based  on  support  diameter  is  of  the  order  of  10  and 
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a  laminar  boundary  layer  will  exist  along  the  length  of  the  support.  If 
negligible  heat  is  lost  from  the  hot-wire  support  itself,  then  the  net 
integrated  heat  transfer  must  be  zero.  This  means  that  any  heat 
transferred  to  the  support  near  the  stagnation  point  will  rapidly  be 
conducted  away  and  returned  to  the  flow  somewhere  downstream.  By 
considering  the  expansion  process  around  a  blunt  body  in  hypersonic 
flow  and  the  size  of  the  stagnation  point  surface  relative  to  the  downstream 
exposed  surface,  it  is  easily  seen  that  the  support  temperature  should  be 
very  close  to  the  adiabatic  temperature  of  a  laminar  boundary  layer  with 
local  external  Mach  number  Mn  (see  the  insert  of  Figure  15).  Expressed 
in  mathematical  terms. 


<T/ro)  % 


1  +  fPT  ( 


tf-i  \  x,  2 

-  >  Mr> 


)  M 


n 


where  M  is  of  order  Z.  5  to  3  for  a  blunt  ogival  body  in  hypersonic  flow. 


n 


This  result  agrees  with  the  experimental  value  to  high  precision. 
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APPENDIX  C 

THE  EFFECT  OF  OVERHEAT  ON  NUSSELT  NUMBER 

During  the  course  of  this  investigation,  a  few  preliminary 

measurements  were  made  to  determine  the  effect  of  overheat  on  the 

30 

measured  Nusselt  number.  From  the  results  of  McClellan  and 

19 

Spangenberg  ,  it  was  not  expected  that  at  low  overheats  and  Mach 
numbers  outside  the  transonic  range  any  significant  effects  would  be 
found.  However,  in  the  region  where  slip  effects  become  appreciable, 
it  is  not  evident  a  priori  that  a  variation  of  surface  temperature  may 
be  neglected.  We  shall  begin  by  discussing  the  data  obtained,  and 
conclude  with  a  qualitative  description  of  the  results  of  this  and  several 
other  investigations. 

Tw  "  Taw 

The  present  data  for  TT  =  - — -  >  0  were  obtained  by  a 

aw 

procedure  identical  to  that  described  in  Section  II.  For  finite  overheat, 
the  variation  of  wire  resistance  with  heating  current  must  be  considered, 
and  Equation  (2)  is  the  appropriate  computational  relation.  Four  tests 
were  carried  out  for  Reynolds  numbers  ReQ  between  0.  5  and  2.  0.  As 
can  be  seen  from  Figures  8,  9,  and  10,  the  effects  (if  any)  of  rarefaction 
should  be  very  noticeable  in  this  region. 

The  data  of  all  four  tests  fell  within  the  shaded  region  shown  in 
Figure  16.  Although  the  reproducibility  of  these  measurements  was 
not  sufficient  to  discern  the  influence  of  Re^  on  the  [Nuq( 't')/NUQ  {Z=  0)J 
relation,  a  comparison  with  the  results  of  McClellan  show  that  even  in 
the  slip  flow  regime,  the  effect  cannot  be  large  for  Mach  numbers  above 
the  transonic  range.  The  difference  between  the  Re 

o 


~  4  data  of 


49 


McClellan  and  the  present  measurements  is  probably  due  to  the  limited 
amount  of  information  available  from  both  sources  (See  McClellan's 
thesis,  p.  19.  ) 

Comparison  of  the  data  of  Spangenberg  at  subsonic  and  transonic 

speeds  with  the  supersonic  and  hypersonic  data  is  facilitated  by  the 

insert  shown  in  Figure  16.  This  insert  (Figure  22a,  Reference  19)  is 

typical  of  the  data  available  in  this  range.  That  the  relation  between 

Nusselt  number  and  overheat  is  quite  complex  may  be  seen  by  comparing 

the  effects  of  7T  for  various  values  of  (M/Re).  (See  for  example 

Figures  22a-/  ,  Reference  19.  ) 

One  fact  which  further  complicates  the  question  of  the  effect  of 

overheat  on  Nuq  is  that  none  of  the  available  experimental  evidence  in 

7  19  21  31 

the  transonic  range  ’  *  ’  includes  the  appropriate  data  for 

zero  overheat.  For  example,  the  data  of  Spangenberg  reproduced  in 

Figure  16  is  normalized  about  the  value  7T  =  0.  214.  As  a  comparison, 
21 

Baldwin  often  found  as  much  as  a  20  per  cent  change  in  the  heat 

transfer  coefficient  between  the  values  "2T  =  0.  07  and  T”  =  0.  22.  Another 

21  7 

rather  surprising  comparison  between  the  works  of  Baldwin  ,  Cybulski  , 
and  Spangenberg  is  the  behavior  found  at  overheats  approaching  one.  As 
illustrated  in  Figure  16,  Spangenberg  found  that  the  change  in  Nu^  was 
close  to  linear  with  a  change  in  overheat  over  a  large  range  of  Mach 
and  Reynolds  numbers.  On  the  contrary.  References  7  and  21  show  a 
negligible  change  in  Nuq  when  1?  is  increased  beyond  about  0.  35. 

In  attempting  to  interpret  these  results,  it  is  necessary  to 
distinguish-between  three  separate  types  of  phenomena  which  are  in 
evidence  in  these  various  investigations.  The  first  may  be  broadly 
classed  as  thermodynamic,  which  would  include  small  changes  in  the 
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local  temperature  of  the  fluid  near  the  body  and  the  subsequent  effects 
on  viscosity,  density,  Prandtl  number  and  the  like.  These  differences 
for  the  most  part  are  small  (at  least  for  small  overheats  or  changes  in 
stagnation  temperature)  and  may  be  accounted  for  at  high  Reynolds 
numbers  by  using  any  of  the  several  available  compressible  boundary 
layer  analyses.  The  second  class  may  be  broadly  termed  induced  flow 
effects,  wherein  a  change  in  body  temperature  is  sufficient  to  alter  the 
basic  features  of  the  flow  field,  and  through  this  "interaction"  change  the 
measured  heat  transfer.  Thirdly  ,  the  effect  of  surface  temperature  on 
molecular  processes  must  also  be  considered. 

From  the  above  delineation,  we  may  immediately  imply  several 
important  generalizations.  First,  in  the  range  of  high  Reynolds  numbers, 
the  effects  of  overheat  should  be  of  the  thermodynamic  type.  Second,  as 

the  Reynolds  number  decreases,  a  given  fractional  change  in  the  boundary 
layer  thickness  caused  by  fluid  property  variation  would  have  an  increasing 

effect  on  the  flow  field  ajacent  to  the  body.  From  boundary  layer  theory, 

_  j, 

q  ~  (d/2)  2  where  d  is  the  diameter  of  the  cylinder.  If  the 

boundary  layer  displacement  thickness  is  not  negligibly  small,  then 

1 


q~  [W2)<1  +  ^>]  ,  and  an  increase  in  the  displacement  thickness 

(i.  e.  ,  an  increase  in  the  wire  temperature)  would  tend  to  decrease  the 

heat  transfer  rate  to  the  body.  This  result  agrees  with  the  experimental 

results  given  in  Figure  16.  The  third  conclusion  is  that  this  induced 

flow  effect  would  be  most  pronounced  in  the  transonic  regime,  where  the 

flow  field  is  extremely  sensitive  to  the  effective  body  shape.  Ail  of  these 

statements  are  borne  out  by  the  present  experiments  and  the  data  of 
19  31 

Spangenberg  and  Winovich  and  Stine  .  In  cases  where  these  different 
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effects  are  of  the  same  order,  it  becomes  difficult  to 

separate  them  by  less  than  rigorous  analysis. 

2  3  6  7  8  Z  Z 

Several  investigations  ’  '  '  ’  *  have  been  conducted  under 
conditions  approaching  free  molecule  flow.  Air,  nitrogen,  helium,  and 
argon  have  all  been  used  as  test  gases  for  these  subsonic  and  supersonic 
experiments.  In  all  cases,  the  test  body  (a  hot-wire,  unheated  wire  or 
butt- welded  thermocouple)  was  maintained  at  a  temperature  differing 
significantly  from  the  equilibrium  recovery  temperature.  Several 
interesting  but  as  yet  unexplained  results  were  obtained  for  various 
combinations  of  overheat  (both  positive  and  negative),  test  medium, 
and  Mach  number.  The  most  important  requirement  at  this  time  is  a 
working  model  of  the  actual  interaction  process  between  a  rapidly 
moving  molecule  and  a  surface  of  arbitrary  roughness  and  temperature. 
As  indicated  earlier,  in  free  molecule  flow  the  body  temperature  will 
change  the  heat  transfer  rate  only  through  its  effect  on  the  accommodation 
coefficient. 
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APPENDIX  D 

MEASUREMENTS  AT  REYNOLDS  NUMBERS 
APPROACHING  FREE  MOLECULE  FLOW 


At  the  conclusion  of  this  investigation,  an  attempt  was  made  to 

j  ,  Af |  yy\  } 

extend  the  data  to  lower  Reynolds  numbers  by  using  .  00005"  diameter  '■ f  ' 

Pt  -  10  per  cent  RH  wires.  These  tests  were  not  quantitatively 
repeatable  because  of  possible  non- uniformity  of  the  wires,  repeated 
breakage  and  consequent  incomplete  calibration,  and  other  reasons 
which  are  not  yet  fully  understood.  However,  each  individual  wire  , 
produced  a  self-consistent  set  of  heat  transfer  data. 

The  accuracy  of  tnese  measurements  was  considerably  less  than 
the  accuracy  obtained  using  .  0001"  diameter  wire.  Microscopic  ex¬ 
amination  revealed  that  these  small  wires  possessed  a  variety  of  small 
"kinks";  whether  these  non- uniformities  were  caused  by  installation, 
tunnel  exposure,  or  the  actual  manufacturing  process  was  not  clear. 

Figure  17  shows  that  the  slope  of  the  Nuq  -  ReQ  curve  obtained 

with  each  individual  wire  was  equal  to  unity  within  the  experimental 

scatter.  Although  the  slopes  agreed  quite  well,  the  proportionality 

constant  varied  by  as  much  as  40  per  cent.  Data  obtained  at  high  Mach 

3 

numbers  in  air  and  nitrogen  by  Weltmann  and  Kuhns  and  by  Wong  are 
included  in  the  figure  for  comparison.  In  each  of  these  experiments, 
the  slope  of  the  Nuq  -  ReQ  curve  is  again  nearly  unity  within  the  limits 
of  the  data.  This  fact  emphasizes  several  of  the  remarks  made  earlier 
with  reference  to  Figure  9;  between  several  independent  investigations 
there  appears  to  be  better  agreement  in  the  slope  of  the  Nuq  -  ReQ 
relation  than  in  the  absolute  magnitude  of  the  correlation. 
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FIG.  4  TYPICAL  HOT-WIRE  INSTALLATION  AFTER  TUNNEL  EXPOSURE 
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FIG.  6  ELECTRICAL  CIRCUIT  DIAGRAM 


Nusselt  Number,  Nu0 


6.0 


4.0 


2.0 

1.0 

0.6 

0.4 

0.2 

0.1 

0.  06 
0.04 

0.02 


Free  molecule  solution 
Stolder, Goodwin  and 

Creoger 

(Pr0  =  0.69,a  =  I.O,M»l) 


0.2  0.4  0.6  1.0  2.0  4.0  6.0  10  20  40  60  100 


Reynolds  Number  ,  Re0 


FIG.  7  EMPIRICAL  CORRELATION  OF  HOT-WIRE  HEAT  TRANSFER  AT 
LOW  REYNOLDS  NUMBERS 


vO 


NUSSELT  NUMBER  —  REYNOLDS  NUMBER  RELATION  FOR  SUPERSONIC 
AND  HYPERSONIC  FLOW 


Slope  of  Nu0  —  Re0  relation, 


very  rotio 


Kn_>  10 


Theoreticol  limit-free-molecule  flow 


o 

Present  Investigotion ,  M  =  5.80 

'M2 

Loufer  ond  Me Clellon, M=  1.33- 4.54 
(without  end  corrections) 

• 

Loufer  and  Me  Glellan,  M  •=  3.05 
(with  end  corrections) 

~  + 

Shermon,  M  =  2.0  and  4.0 

X 

Stalder, Goodwin  and  Creager 

[IjLjl 


E  xperimentol  limit  —  continuum  flow 


D.O  2 

0.04  0.06 

0.10 

0.2 

0.4  0.6  1.0 

2.0  4.0  6.0 

Free-streom  Knudsen  number, 

Knco 

FIG.  10 

NORMALIZED 

VARIATION 

OF 

RECOVERY  TEMPERATURE 

WITH  KNUDSEN  NUMBER 

Wire  resistance,  ohms 


FI6.II  -DATA  OBTAINED  FROM  CYLINDER  WAKE  TRAVERSE 


oo 


c)  Velocity  defect 


Correction  quantity  , 


Measured  support  temperature  , 


0  5  10  20  30  40  60  80  100 

Stagnation  pressure,  p s i  g 


FIG.  15  WIRE  SUPPORT  TEMPERATURE  AS  A  FUNCTION  OF  TUNNEL  PRESSURE 
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Attention:  Aeroballistics  Division,  M-AERO-E,  Mr.  T.  G.  Reed  (3  copies) 

NASA 

Langley  Research  Center 
Langley  Field,  Virginia 
Attention:  Librarian 

Attention:  Mr.  Clinton  E.  Brown,  Chief,  Theoretical  Mechanics  Division,  Bldg.  121 
Attention:  Dr.  Adolf  Busemann 

Attention:  Mr.  Cnarles  H.  McLellan,  1 1 -Inch  Hypersonic  Tunnel  Section 
NASA 

Ames  Research  Center 
Moffett  Field,  California 
Attention:  Library 

NASA 

Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland  35,  Gnio 

Attention:  Library,  Mr.  George  Mandel 
2  copies 

Miscellaneous  Government  Agencies 

United  States  Atomic  Energy  Commission 
P.  O.  Box  62 
Oak  Ridge,  Tennessee 
Attention:  Library 

U.  S.  Department  of  Commerce 
National  Bureau  of  Standards 
Washington  25,  D.  C. 

Attention:  Dr.  G.  B.  Schubauer 

Chief,  Fluid  Mechanics  Section 


Universities  and  Non-Profit  Organizations 


Brown  University- 
Division  of  Applied  Mathematics 
Providence  12,  Rhode  Island 
Attention:  Professor  R.  E.  Meyer 

Brown  University 
Division  of  Engineering 
Providence  12,  Rhode  Island 
Attention:  Dr.  Ronald  F.  Probstein 

University  of  California  at  Berkeley 
Aeronautical  Sciences  Department 
Room  203,  Mechanics  Building 
Berkeley  4,  California 
Attention:  Professor  S.  A.  Schaaf 

University  of  California 
405  Hilgard  Avenue 
Los  Angeles  24,  California 
Attention:  Engineering  and  Mathematical 
Sciences  Library 

University  of  California 
Department  of  Engineering 
Los  Angeles  24,  California 
Attention:  Professor  A.  F.  Charwat 
Attention:  Professor  John  W.  Miles 
Attention:  Dr.  N.  Rott 

Case  Institute  of  Technology 
Department  of  Mechanical  Engineering 
University  Circle 
Cleveland  6,  Ohio 
Attention:  Dr.  G.  Kuerti 

Columbia  University 

Department  of  Mechanical  Engineering 
New  York  27,  N.  Y. 

Attention:  Professor  Robert  A.  Gross 
Cornell  University 

Graduate  School  of  Aeronautical  Eng. 
Ithaca,  New  York 
Attention:  Library 
Attention:  Dr.  William  R.  Sears 

University  of  Florida 
Department  of  Aeronautical  Eng. 
Gainesville,  Florida 

Attention:  Professor  David  T.  Williams 

Harvard  University 
Division  of  Eng.  and  App.  Physics 
Cambridge  38,  Massachusetts 
Attention:  Dr.  Howard  W.  Emmons 


University  of  Illinois 

Department  of  Aeronautical  Engineering 
Urbana,  Illinois 

Attention:  Professor  Harold  O.  Barthel 
Attention:  Dr.  Allen  I.  Ormsbee 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attention:  Dr.  L.  L.  Cronvich 
Attention:  Dr.  F.  K.  Hill 

The  Johns  Hopkins  University 
Department  of  Mechanics 
Baltimore  18,  Maryland 
Attention:  Dr.  Francis  H.  Clauser 
Attention:  Dr.  Stanley  Corrsin 
Attention:  Professor  L.  S.  G.  Kovasznay 

Lehigh  University 
Department  of  Physics 
Bethlehem,  Pennsylvania 
Attention:  Dr.  Raymond  J.  Emrich 

University  of  Maryland 

Department  of  Aeronautical  Engineering 
College  Park,  Maryland 
Attention:  Professor  S.  F.  Shen 

University  of  Maryland 
Institute  for  Fluid  Dynamics  and 
Applied  Mathematics 
College  Park,  Maryland 
Attention:  Director 
Attention:  Professor  J.  M.  Burgers 
Attention:  Professor  Francis  R.  Hama 
Attention:  Professor  S.  I.  Pai 

Massachusetts  Institute  of  Technology 
Department  of  Aeronautics  and  Astronautic 
Aerophysics  Laboratory 
560  Memorial  Drive 
Cambridge  39,  Massachusetts 
Attention:  Dr.  Morton  Finston 

Massachusetts  Institute  of  Technology 
Department  of  Mathematics 
Cambridge  39,  Massachusetts 
Attention:  Professor  C.  C.  Lin 
Attention:  Dr.  George  B.  Whitham 


Massachusetts  Institute  of  Technology 
Department  of  Mechanical  Engineering 
Cambridge  39,  Massachusetts 
Attention:  Dr.  A.  H.  Shapiro 
Attention:  Dr.  H.  Guyford  Stever 
Room  3-174 

Massachusetts  Institute  of  Technology 
Department  of  Aeronautics  and  Astronautics 
Cambridge  39,  Massachusetts 
Attention:  Professor  E.  Mollo-Ghristensen 
Room  33-412 

Attention:  Dr.  Leon  Trilling 
Room  33-412 

University  of  Michigan 
Ann  Arbor,  Michigan 
Attention:  Engineering  Library 

University  of  Michigan 
Willow  Run  Laboratories 
P.  O.  Box  618 
Ann  Arbor,  Michigan 
Attention:  BAMIRAC  Library 

Mr.  Richard  Jamron,  Head 
Information  Handling  Group 

University  of  Michigan 

Aeronautical  Engineering  Laboratories 

Aerodynamics  Laboratory 

North  Campus 

Ann  Arbor,  Michigan 

Attention:  Mr.  James  L.  Amick 

University  of  Michigan 

Aeronautical  Engineering  Laboratories 

Aircraft  Propulsion  Laboratory 

North  Campus 

Ann  Arbor,  Michigan 

Attention:  Professor  J.  A.  Nicholls 

University  of  Michigan 
Department  of  Aeronautical  and 
Astronautical  Engineering 
Ann  Arbor,  Michigan 
Attention:  Dr.  Arnold  M.  Kuethe 
Attention:  Professor  V.  C.  Liu 
Attention:  Piofessor  William  W.  Willmarth 

University  of  Michigan 
Department  of  Physics 
Ann  Arbor,  Michigan 
Attention:  Dr.  O.  Laporte 


University  of  Minnesota 

Institute  of  Technology 

Rosemount  Aeronautical  Laboratories 

Rosemount,  Minnesota 

Attention:  Mrs,  Lmda  Caldon, 

Librarian 

New  York  University 

Institute  of  Mathematics  and  Mechanics 
53  Washington  Square,  South 
New  York  12,  N.  Y. 

Attention:  Library 

North  Carolina  State  College 
Department  of  Mechanical  Engineering 
Raleigh,  North  Carolina 
Attention:  Professor  R.  M.  Pinkerton 

Northwestern  University 
The  Technological  Institute 
Evanston,  Illinois 

Attention:  Professor  Ali  Butent  Cambel 

The  Ohio  State  University 
Department  of  Aeronautical  and 
Astronautical  Engineering 
2036  Neil  Avenue 
Columbus  .10,  Ohio 
Attention:  Professor  John  D.  Lee 
Attention:  Professor  Gavin  L.  VonEschen 

Polytechnic  Institute  of  Brooklyn 
Aerodynamics  Laboratory 
527  Atlantic  Avenue 
Freeport,  New  York 
Attention:  Library 

Attention:  Professor  Martin  H.  Bloom 
Attention:  Professor  Antonio  Ferri 
Attention:  Professor  Paul  A.  Libby 

Princeton  University 
School  of  Engineering 
James  Forrestal  Research  Center 
Princeton,  New  Jersey 
Attention:  Library 

Attention:  Gas  Dynamics  Laboratory 
Attention:  Dr.  Seymour  Bogdonoff 
Attention:  Professor  Sin-I  Cheng 
Attention:  Dr.  Luigi  Crocco 


Purdue  University 
School  of  Aeronautical  and 
Engineering  Sciences 
Wast  Lafayette,  Indiana 
Attention:  Aero,  and  Engineering 
Sciences  Library 

Rensselaer  Polytechnic  Institute 
Department  of  Aeronautical  Engeering 
Troy,  New  York 
Attention;  Library 
Attention:  Dr.  Ting-Yi  Li 

University  of  Rochester 

College  of  Engineering 

Department  of  Mechanical  Engineering 

River  Campus  Station 

Rochester  20,  New  York 

Attention:  Professor  Martin  Lessen 

University  of  Southern  California 
Engineering  Center 
University  Park 
Los  Angeles  7,  California 
Attention:  Director 
Attention:  Dr.  H.  T.  Yang 

University  of  Southern  California 
Engineering  Center 

Aeronautical  Laboratories  Department 
P.  O.  Box  1001 
Oxnard,  California 
Attention:  Mr.  J.  H.  Carrington, 
USCEC-ATL 

Stanford  University 

Department  of  Aeronautical  Engineering 
Stanford,  California 

Attention:  Professor  Daniel  Ber shader 
Attention:  Dr.  Milton  Van  Dyke 
Attention:  Professor  Walter  G.  Vincenti 

University  of  Texas 

Defense  Research  Laboratory 

P.  O.  Box  8029 

Austin  12,  Texas 

Attention:  Dr.  M.  J.  Thompson 


University  of  Washington 

Department  of  Aeronautical  Engineering 

Seattle  5,  Washington 

Attention:  Engineering  Librarian 

Attention:  Professor  R.  E.  Street 

University  of  Wisconsin 

Theoretical  Chemistry  Laboratory 

P.  O.  Box  2127 

Madison  5,  Wisconsin 

Attention:  Dr.  Joseph  O.  Hirschfelder 

Yale  University 

Department  of  Mechanical  Engineering 
New  Haven,  Connecticut 
Attention:  Dr.  Peter  Wegener 
Attention:  Dr.  Alan  Kistler 

Institute  of  the  Aerospace  Sciences 
2  East  64th  Street 
New  York  21,  New  York 
Attention:  Library 


University  of  Virginia 
Department  of  Physics 
Charlottesville,  Virginia 
Attention:  Dr.  Jesse  W.  Beams 


Industrial  and  Research  Companies 


Aeronautical  Research  Associates 
of  Pridceton,  Inc, 

50  Washington  Road 
Princeton,  New  Jersey 

Attention:  Dr.  Coleman  dup.  Donaldson 
Aeronutronic 

A  Division  of  Ford  Motor  Company 

Ford  Road 

P.  O.  Box  697 

Newport  Beach,  California 

Attention:  Dr  L.  L.  Kavanau 

Advanced  Programs  Staff 


Aerospace  Corporation,  Inc. 

P.  O.  Box  95085 
Los  Angeles  45,  California 
Attention:  Dr.  Chieh-Chien  Chang 
Attention;  Dr.  J.  Logan,  Director 
2  copies 

Attention:  Dr.  H.  Mirels 


ARO,  Inc. 

Arnold  Air  Force  Station 
Tennessee 

Attention:  AEDC  Library 
Attention:  Dr.  B.  H.  Goethert 

Director  of  Engineering 
Attention:  TS(T1) 

ARO,  Inc. 

von  Karman  Gas  Dynamics  Facility 
Arnold  Air  Force  Station 
Tennessee 

Attention:  Dr.  J.  Lukasiewicz,  Chief 
Attention:  Mr.  J.  Leith  Potter,  \ 

Manager,  Research  Branch 

AVCO-Everett  Research  Laboratory 
2385  Revere  Beach  Parkway 
Everett  49,  Massachusetts 
Attention:  Barbara  A.  Spence, 

Technical  Librarian 

AVCO  Research  and  Advanced 
Development  Division 
201  Lowell  Street 
Wilmington,  Massachusetts 
Attention:  Mr.  A.  Kahane 

Assistant  Technical  Director 
Attention:  Dr.  Frederick  R.  Riddell, 

Tech.  Ass't.  to  Pres.  -S.  Tec. 


Boeing  Airplane  Company 
Aero-Space  Division 
Seattle  24,  Washington 
Attention:  Library  13-84 

CONVAIR 

A  Division  of  General  Dynamics  Corp. 

Astronautics  Division 

P.  O.  Box  1128 

San  Diego  12,  California 

Attention:  Mr.  K.  J.  Bossart,  Tech.  Dii 

Attention:  Mr.  W.  B.  Mitchell,  595-10 

CONVAIR 

A  Division  of  General  Dynamics  Corp. 
Scientific  Research  Laboratory 
5001  Kearny  Villa  Road 
San  Diego  11,  California 
Attention:  Mr.  William  H.  Dorrance 
Senior  Staff  Scientist 
Mail  Zone  1-162 
Attention:  Mr,  Merwin  Sibulkin 
Staff  Scientist 

CONVAIR 

A  Division  of  General  Dynamics  Corp. 
San  Diego  12,  California 
Attention:  Dr.  H.  Yoshihara 

Chief  of  Fluid  Dynamics  Res. 
Mail  Zone  6-105 

CONVAIR 

A  Division  of  General  Dynamics  Corp. 
Aerospace  Technology  Section 
P.  O.  Box  748 
Fort  Worth,  Texas 
Attention:  Mr.  R.  C.  Frost 
Dept.  6-1 
Mail  Zone  E63 

CONVAIR 

A  Division  of  General  Dynamics  Corp, 
Fort  Worth,  Texas 
Attention:  Mr.  A.  P.  Madsen 

Aerodynamics  Group  Eng. 
Mail  Zone  E63 

Attention:  Mr.  W.  G.  McMullen 
Attention:  Mr.  Robert  H.  Widmer 

CONVAIR 

A  Division  of  General  Dynamics  Corp. 

Daingerfield,  Texas 

Attention:  Mr.  J.  E.  McMichael 

Chief,  Jet  Engine  Department 


Cornell  Aeronautical  Laboratory,  Inc. 

P.  O.  Box  235 
Buffalo  21,  New  York 
Attention:  Library 
Attention:  Dr.  A.  H.  Flax 
Attention:  Mr.  A.  Hertzberg 

Head,  Aerodynamic  Research 

Douglas  Aircraft  Company,  Inc. 

Missiles  and  Space  Systems 
3000  Ocean  Park  Blvd. 

Santa  Monica,  California 
Attention:  Library 
Chief, 

Aero/Astrodynamic s  Section 
2  copies 

Attention:  Mr.  R.  J.  Gunkel 
Chief, 

Aero/Astrodynamic s  Section 

Douglas  Aircraft  Company,  Inc. 

827  Lapham  Street 
El  Segundo,  California 
Attention:  Dr.  A.  M.  O.  Smith 


Lockheed  Aircraft  Corporation 
Lockheed  Missiles  and  Space  Company 
Technical  Information  Center  (50-14) 

3251  Hanover  Street 
Palo  Alto,  California 
Attention:  Dr.  W.  A.  Kozumplik 
3  copies 

Lockheed  Aircraft  Corporation 
Lockheed  Missiles  and  Space  Company 
P.  O.  Box  504 
Sunnyvale,  California 
Attention:  Mr.  R.  Smelt, 

Chief  Scientist 

Lockheed  Aircraft  Corporation 
Lockheed  Missile  and  Space  Company 
Palo  Alto,  California 
Attention:  Mr.  Maurice  Tucker 

Spacecraft  and  Missiles  Research 


General  Electric  Company  (MSVD) 
Space  Technology  Center 
Space  Sciences  Laboratory 
King  of  Prussia,  Pennsylvania 
Attention:  Dr.  H.  Lew 


General  Electric  Company 
Research  Laboratory 
P.  O.  Box  1088 
Schenectady,  New  York 
Attention:  Dr.  Henry  T.  Nagamatsu 

Giannini  Controls  Corporation 
1600  South  Mountain  Avenue 
Duarte,  California 
Attention:  Library 

Grumman  Aircraft  Engineering  Corp. 

Bethpage,  New  York 

Attention:  Mr.  Charles  Tilgner,  Jr. 

Hughes  Aircraft  Company 
Culver  City,  California 
Attention:  Mr.  E.  O.  Marriott 

Manager,  Aerodynamics  ] 


Lockheed  Aircraft  Corporation 
Missiles  and  Space  Division 
7701  Woodley  Avenue 
Van  Nuys,  California 
Attention:  Library 

Lockheed  Aircraft  Corporation 
Marietta,  Georgia 
Attention:  Dr.  W.  F.  Jacobs 

Aerodynamics  Dept.  -72-07 

Marquhardt  Aircraft  Company 
P.  O.  Box  2013  -  South  Annex 
Van  Nuys,  California 
Attention:  Technical  Library 

The  Martin  Company 
Baltimore  3,  Maryland 
Attention:  Mr.  K.  Jarmolow 
Mail  No.  J-3033 

Attention:  Dr.  Mark  V.  Morkovin 
Mail  No.  J-3033 

McDonnell  Aircraft  Corporation 
Lambert- Saint  Louis  Municipal  Airport 
P.  O.  Box  516 
St.  Louis  66,  Missouri 
Attention:  Mr.  Kendall  Perkins 


The  RAND  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Attention:  Librarian 

Attention:  Dr.  Carl  Gazley,  Jr. 

Attention:  Mr.  E.  P.  Williams 

Aero-Astronautics  Dept. 


Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  New  York 
Attention:  Engineering  Library 
Attention:  Mr.  R.  W.  Perry 

Chief,  Re-Entry  Simulation  Lab. 
Applied  Research  and  Dev. 


Space  Technology  Laboratories,  Inc. 
P.  O.  Box  95001 


Los  Angeles  45,  California 
Attention:  Technical  Information  Center 
Document  Procurement 
Bldg.  C,  Room  2412 
Attention:  Dr.  James  E.  Broadwell 

Aerodynamics  Research  Section 
Attention:  Dr.  C.  B.  Cohen 
Attention:  Dr.  Luois  G.  Dunn,  President 
Attention:  Dr.  Andrew  G.  Hammitt,  Head 
Aerodynamics  Research  Section 
Attention:  Mr.  Ernest  I.  Pritchard 
Attention:  Dr.  George  E.  Solomon 


Internal 


Mr.  Paul  E.  Baloga 
Dr.  Julian  D.  Cole 
Dr.  Donald  E.  Coles 
Dr.  Anthony  Demetriades 
Dr.  Toshi  Kubota 
Professor  Lester  Lees 
Dr.  H.  W.  Liepmann 
Dr.  Clark  B.  Millikan 
Dr.  Barry  L.  Reeves 
Dr.  Anatol  Roshko 

Dr.  Frank  Marble 
Dr.  S.  S.  Penner 
Dr.  W.  D.  Rannie 
Dr.  Edward  Zukoski 

Dr.  Harry  Ashkenas 
Mr.  George  Goranson 
Dr.  James  M.  Kendall 
Dr.  John  Laufer 
Dr.  Thomas  Vrebalovich 
Mr.  Richard  Wood 

Aeronautics  Library  (2) 

Hypersonic  Files  (3) 

Hypersonic  Staff  and  Research  Workers  (20 


Sperry  Utah  Engineering  Laboratory 
Division  of  Sperry  Rand  Corporation 
322  North  21st  Street  West 
Salt  Lake  City  16,  Utah 
Attention:  Mr.  Malcolm  L.  Matthews 


Systems  Corporation  of  America 
1007  Broxton  Avenue 
Los  Angeles  24,  California 
Attention:  Dr.  Paul  D.  Arthur 

United  Aircraft  Corporation 
Research  Laboratories 
East  Hartford  8,  Connecticut 
Attention:  Mr.  John  G.  Lee 


